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Chapter 1. Eigenvalues and Eigenvectors for Eybrld Coordinate 
Equations of Motion for Flexible Spacecraft 

ABSTRACT . Literal characterizations are developed for the eigenvalues and 
eigenvectors of a system of linear time-invariant equations which describes 
the attitude motion of flexible spacecraft In terms of hybrid coordinates. 

The elgenproblem Is shown to reduce to that of a symmetric, and positive 
definite matrix of lower dimension. Both analytical and mlnlrnax characteri- 
zation methods prove to be useful In localizing the eigenvalues for the zero 
damping case. A perturbation method Is employed to Investigate th.e effects of 
modal damping. The resulting eigenvectors generate a canonical form, based 
on which the controllability is examined. 
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1. INTRODUCTION 


One of the most commonly accepted procedures for the description of the 
attitude dynamics and control of flexible spacecraft Is a hybrid coordinate 
formulation, In which the attitude variables of a primary body are discrete 
coordinates while the displacements of elastic appendages relative to the 
primary body are characterized by distributed or modal coordinates. Since 
the formulation has been established and the linearized variational equations 
are at hand,[l] It concerns us greatly to examine the system behavior based 
on this representation. In this examination the elgenproblems play a key role 
for many reasons. Firstly, the eigenvalues or eigenvectors represent essential 
features of the system such as stability, steady state response and so on. 
Moreover, in application of existing control theory, whether It be modern 
control theory or classical, eigenvalues and eigenvectors are useful In 
canonical transformation, pole location and so on. Secondly, knowledge of 
elgenpropertles of the system affords a sound basis for the truncation 
procedure which Is essential to system simulation or control system design. 

However, such systems are usually multivariable systems with an extremely 
high dimension, so that numerical calculation of the eigenvalues or eigen- 
vectors Is not an easy task. Hence, It Is desirable to characterize them In 
some closed form. 

The purpose of this paper Is to characterize the eigenvalues and eigen- 
vectors of the system in literal expressions, by utilizing peculiar properties 
of the system parameter matrices as they have been previously established. 

The eigenvalues are localized in terms of inertial matrices and modal param- 
eters, and a procedure for calculation of the eigenvectors is proposed in some 
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special cases of practical Interest. The effects of truncation and damping 
are also examined. 

A canonical form of the system state equations Is derived, and 
controllability criteria are discussed and compared with the previous 
results. [3] 

In the hybrid coordinate representation, the vehicle translational 
equations take the form[l] 

1*0 - 6% = T (la) 

and the appendage deformation equations" may be written [1] ,[ 3] 

•• • 2 •* T 

n + 2can + o n - 60 = i? t (lb) 

where T Is the 3x1 matrix representing the external torque vector about the 
vehicle mass center c, for an orthonormal vector basis in the primary body b, 
I* is the 3x3 inertia matrix of the total vehicle for c, 0 is the 3x1 
matrix of 1-2-3 Inertial attitude angles of b, and establishes the location 
and type of the attitude control actuators. [3] Here r) is the N x 1 matrix 


representing the modal coordinates, where N represents the number of appendage 
modes; (|) is an N x 6n matrix whose columns <t>^ represent the mode shapes of 

/ 

those appendage vibrations with nonzero natural frequencies , which could 
occur independently of each other If the primary body were translationally free 
but constrained against rotation; O is an N x N diagonal matrix whose elements 


are arranged in non-decreasing order; and ^ is also an N x N diagonal 
matrix whose elements represent the modal damping associated with modal 
frequencies for j = 1, ...» N. The N x 3 matrix 6 is established by ({) 
and the geometry and mass distribution characteristics of the appendages, as 
in Eq. (278) of [1]. It should be noted that (j> is normalized with respect to 
the generalized inertia matrix, as in Eq. (213) of [1], and that with this 
normalization 6 has the dimension of square root of inertia. 
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For convenience In later discussions, we define the following matrices: 
Y A (I*)*®0 (2j 

A A 6 (I*)”*® (21 

u A (I*)"S: (2< 


D A 2t,a 
8 

With the definitions of Eqs. C2) , Eqs, Cl) become 

•• _.«• 

T 

Y - A ri = u 


ri + Dn + cj^n - Ay = (i*)**’u 


where 


with s 02 5 . . . sOjj. 

Furthermore, it is convenient to decompose A into submatrices 


where the 1x3 matrix A'^ is defined by 

A^ A [ A^^, A^2* ] . j = 1, 2, N 

and the scalar quantities A^j^ (j = 1,2,..., N; k = 1,2,3) 
are dimensionless and [1] 

0 < A^, <1 


t Since I* is symmetric and positive definite by nature, (I*)^ is always 

feasible with the following interpretation: there exists an orthogonal 

T D D 

matrix R such that RI*R = I where I is a diagonal matrix with positive 
elements, and (^I)^ is defined. This permits I* = R^ ^IR = R^(^I)^ RR^(^I)Sl, 
so that (I*)^ = R^(°I)Sl 
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T 

In what follows, the 3 x 3 matrix A A plays an Important role so that some of 
its properties are stated; if we interpret 6 in terms of primitive definitions, 
then we obtain the physical interpretation.[ 2 ] 

= I* - 1° 

0 

where I is the Inertia matrix of the primary body referred to its own mass 
center. From Eq. (2b), 

a'^A = (Ad) 

0 T 

Obviously, the matrix (I*-I ) Is: symmetric and positive definite and hence A A 

is positive definite. Furthermore, if we consider the matrix 

■ . : 1 ! '■ ' t . ( t 5 1 ' i . - ■ i 

U 3 - A^A = (I*)"’®I°(I*)"^ (4e) 

T 

it Is also positive definite. Therefore, the eigenvalues of A A are greater 
than zero and less than unity, provided that there exist at least three 
independent rows in A (see Appendix A) . 

In terms of the (2N + 6 ) = x 1 state variable 



Eqs. (3) may be written as 
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with 


and 


82 = M2 + (I*)^ 

B. 6 AM„ + M, (p'^ ^(I*)^ 

4 2 1c 


Ml 4 (U^ - AA^)-^ 


M2 ^ (U3 - a'^a)'^ 


(6c) 

(6d) 

(6e) 

(6f) 


T T 

The noted properties of A A guarantee the nonslngularlty of (U^ - A A) , so 


that M2 Is always feasible. Since |u^ - AA 


I U3 - A^A I , Ml Is also 


1 1 . i ; f i 


1 . 


feasible. Moreover, Mi and M2 are both symnetrlc and positive definite, as 

it. 

shown In Appendix A. 


2. REDUCTION OF EIGENPROBLEM 


Let X be a generic eigenvector of A corresponding to a generic eigenvalue 
X; then the problem Is to find X and x which satisfy 


l*“a(+6 - AI » - » 


C7) 


with the characteristic equation 


^^2N+6 ~ ° 


( 8 ) 


By the determinant partitioning formula [6], Eq. (8) may be written from Eq. 
(6a) as 


I 0 


XU- 

-u. 


XU„ 

-u„ 

3 

3 


N 

N 

0 

XU3 


MiO^ 

XU„ + M,D 
N 1 


= 0 


(9) 


From the first determinant of Eq. C9) 


X^ = 0 


( 10 ) 


and we have X = 0 with multiplicity 6. 
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The second determinant of Eq. (9) is further decomposed Into 

v*“l> 


or 


|XUjj|-|AUjj + M^D - I = 0 


IX^, + XM.D + M.a^l = 0 
N 1 i 


( 11 ) 


If we also decompose the eigenvector x as 
X, 


L*4 J 


( 12 ) 


where Xj^,X 2 ,X 2 and x^ are 3x1, 3x1, Nxl, and Nxl matrices associated respec- 
. rtlvely with Yt1'»n and n» then Eq. (7) becomes 


XUj -?3 

0 XU. 


[;:1 


T 2 T 

A A Mj^D 


[:;] ■ • 


(13) 


and 




-U„ 


XUj^+Mj^D 


L *4 J 


(14) 


This decomposition Indicates a procedure for obtaining the eigenvectors: 
since Eq. (14) does not Include Xj^ or x^> we can first solve Eq. (14) for x^ 
and x^, and then substitute the solution' Into Eq. (13) to produce x^^ and X 2 . 
For the solution X » 0 of Eq. (10) , the determinant of Eq. (14) becomes 

I. . ' - 

= (-1)*^ 0 


0 -Uj, 

. M^D 


and hence there exists only the trivial solution for x^, x^, l.e.. 
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[::] 


( 15 ) 


Substituting Eqs. (10) and (15) into Eq. (13) yields 

0 -U. 

^ I I I 

= 0 


Implying 


= arbitrary, x^ = 0. 


I (1) (2) 

We may introduce three independent and orthogonal 6x1 matrices C , C , and 

1 i 

defined by 


.( 1 ) _ 


* 1 ' 

i - * 

■ 0. ■ 


’ o' 

0 


1 

* 

0 

0 


. . 0, , 


1 

0 


0 

f 

0 

0 


0 

' 1 

0 

. 0 . 


. 0 . 


.0 . 


r :-(l6a) 


The six eigenvectors corresponding to X = 0 are hull except for the upper 6x1 
partitions, which are and three linear combinations of these 

matrices. We see that the number of independent eigenvectors corresponding 

I 

to X=0 is less than the multiplicity 6 of X = 0, so the matrix A is not similar 
to a diagonal matrix but to a nondiagonal Jordan form. We can find such a 
transformation matrix if we construct the generalized eigenvectors, which are 
null except for upper 6x1 partitions given by (among many) 


(16b) 



r d 1 




ro 1 


0 


0 


0 


0 

1 

. 

-0 

0 

and g(3) = 

0 

0 


0 


1 


0 


L 0 J 


L 0 J 


L 1 J 


For the nonzero eigenvalues determined by Eq. (11), we have from Eqs. (13) 
and (14) 

(M.O^ + XM n + X^U„) x, = 0 (17a) 

1 1 N j 
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(17b) 


1 T 2 T 

- y (A + A x^) 


(17c) 


(17d) 


Equations (17) indicate that x^ is the only vector to be solved as the eigen- 
vector, and that other vectors are all determined by matrix or scalar multi- 
plications and additions. 

Thus, the eigenvalue and eigenvector problem of A reduces to Eq. (17a), 
in which A is in the second order form. Eq. (14) is also a reduced form vrfilch 
may be written in a standard form of eigenvalue problems as 


where 


(AU£j, - Y = 0 


A 0 ' 

.V 4 {. _N_ 

. ( 


C19a) 


[:;] 


C19b) 


In order to localize the eigenvalues and characterize the eigenvectors , 
we will take the following procedures: 

We first assume that ^ = 0 (so D = 0) and in the second order form of 
Eq. (17a) we define 
A 2 

y = -A^ (2C 

to which corresponds the eigenvector \p. Then Eq. (17a) becomes 

(yUjj - =0 (21 


for which the noted properties of will be fully utilized. 
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After obtaining the results for D = 0, we will return to the first order 
form of Eq. (lA) and treat D as a small perturbation. 


3. CHARACTERIZATION OF EIGENVALUES 


3.1 Analytical Method 

Since is nonsingular, we may define the characteristic polynomial, 
f(y) , for Eq. (21) by j 


f(y) = M^O^I 

or, in view of Eq. (6e) , 


(22a) 


f(y) = |yUjj - o-^^yAA^I (22b) 

If we rewrite f(y) as 

f(y) = |m‘^| la| |yUjj - M^a^l |o“^| 

= |m"^ 1 |yUjj - 0M^0| 

then w^'recognize that y is the eigenvalue of an NXN real, symmetric, and 
positive definite matrix 0M^^0. Thus, we have[4] the following fact. 

Fact 1 . All the roots of f(y) = 0 are real, and positive, i.e., all the 
2 

eigenvalues of M^^0 are real and positive, and hence, by Eqs. (10) and (20), 

all the nonzero eigenvalues of are imaginary. In addition, the eigenvalue 

problem is well-conditioned (a small change in the elements of 0M^0 does not 

cause any abrupt change in the eigenvalues) . 

In what follows we attempt to characterize the eigenvalues further mainly 

2 

based upon a particular structure of f(y) (Eq. (22b)), in which yU^j - 0 is 
T 

diagonal and AA is of rank 3. 

In order to do so, we first tentatively assume that for the root of 
2 

f(y) = 0, (yU^ - 0 ) is nonsingular, i.e., : 

|yUjj - 0^N 0 (23) 
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Under this assumption, we may rewrite Eq. (22) using the determinant identity [ 6 ] 
as 

f(li) =. lyUjj - 0^1 • g(y) (24a) 

where 

g(M) = IU 3 - (pU^ - a^)"^A| (24b) 

Since in this case the roots of f(jj) = 0 are identical to those of g(y)=0, 

we will treat g(y) rather than f(p) because the former is related to a 35<3 

matrix while the latter is to an NxN matrix. 

The cases in which Eq. (23) is violated will be discussed in such a way 

as to establish under what conditions the system eigenvalue becomes 
2 

identical to for some Z. A sufficient condition is stated below for this 
to hold. 

Z 2 

Fact 2 . If A =0, then = 0. This means that the natural frequency, 

£ 

cr^, of an appendage vibration mode with A = 0 is also that of a system vibra- 
tion mode. 


Proof. If we rewrite Eq. (22b) as 


f(U) = 


IJ 

1 1 X 


aIaNT. 

-y A A - 
Ir 

I 


I .A. IT \ 2 .Z.ZT ''a^aNT 

-y A A y -o» -y A A ^ -y A A 

! ^ . "v. I 

' aNaI 2 I .N.NT 

-yAA y -o„ -y A A 

N 


’ z 

and expand about the. )l-th row with A = 0, we have 


fCy) = Cuyp 


2 a^aIT 
y -yA A ^ 

1 

aNaIT 

-yA A - 


aIaNT 

-yA.A 


N.NT 


•y-Ojj -yAA 
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so that 


f(op = 0 . 


It is noted that = 0 does not necessarily imply a = 0. In fact, if 

1 2 
Oj's and A'^’s satisfy a certain equality, then it happens that = 0 for 

£ 

nonzero A as will be shown later. 

In order to evaluate the roots of Eq. (24b), we define 

A 


dj(ii) 


J = 1,2,. ...N 


\i - a 


j 


(25) 


and, with this definition, we have 


y(yUjj -ch ^ 


dj^(y) 


d2(y). 




(26) 


and 


g(y) = 


^3- 


N 


T ^ 


d,(y) 

^ S 

N 


LA' 


U - E d (m)A^^A^ 

d j=l J 

N 

i 


rii 


d^Cy) J La*^J 


(27) 


C( 

where 6^ is the Kronecker delta defined as 


- 

S - 


1 if a = B 

0 if a B 
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and the quantity In the brackets Indicates the (a, 3) element of 
[U^ - (J A (p - o ) A]. Some properties of dj(y) are given In Appendix B, 
The function g(y) is easily evaluated for some particular values of 
interest. Since dj(y) = 0 for y = 0, 


g(0) = |UJ = 1. 


(28) 


and since 


lim d.(y) = 1. 
y-H:® ^ 


lim g(y) = |U. - A^A| 
y>±«> 


(29) 


Consequently 


0 < lim g(y) < 1 


(30) 


in view of the noted property of Eq. (4e) . 

Considering the properties of g(y) of Eqs. (28) and (30) , we will further 

characterize the roots of g(y) = 0 based on the fact that it grows without 

2 

bound when y approaches to for any j = 1,2,...,N. 

As proven in Appendix C, g(y) is expressed by 

N 


where 


g(y) = 1 - 13 G (y) dj(y) 
S,=l ^ ^ 


G (y) = - E d (y) kJ 

^ j=l 3 3 


(31) 


iSj <kSJl 




(32) 
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with 


and 


where 


kJ = (A^'S^KA^)'^ 


= j A^CA^S:^)"^ j 


(33) 


(34) 


A^ = 


0 -M 




0 -ti 




(35) 


Z i 

We recognize the following properties of K, and L., , noting that by the 

j 

assumption in Eq. C23) and Fact 2 following Eq. C24b) we have excluded th.e 
0 

case A = 0 for all 


(i) kJ 20 (36) 

with the equality when A A^ = 0> Implying that two vectors have 
the same direction; 


(ii) kJ = kJ 


(v) 


(37) 


(iii) 20 

with the equality when A^(A^A^ =0. This takes place if any two 

1k f- 

of A*^, A and A have the same direction or if one of these three 
is perpendicular to the cross-product of the others. Hence 

(iv) If k!" = 0, then lJ = 0. (38) 

SL 


“jk 




Ak 


(39) 


From the expression of Eq. (31) with the definition of Eq. (32) and the 
noted properties of Eqs. (36) -(39) we have the following fact: 
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2 

Fact 3. All the roots of gCy) = 0 are greater than . 

2 

Proof. Consider the value of gCh) for 0 < y < Then, by the definition 

of Eq. (25) 

^ — 2 ^ ° ^ ^ (^Oa) 

and hence 

Gj^(y) i 0 lor £ = 1,2 N (40b) 

Equations (40) guarantee that 

2 

g(y) a 1 for 0 < y < 

2 

so that there exist no roots of g(y) = 0 in the range 0 < y < . On the 

other hand, from Fact 1, there exists no negative root of g(y) = 0. Therefore, 

2 

all the roots of g(y) = 0 are greater than a^. 

One can interpret the Facts 1 and 3 in terms of the eigenvalues, X, of A 
as follows. If the roots of g(y) = 0 are denoted by y^ , j = 1,2,...,N and 
arranged in non-decreasing order, then the Fact 3 indicates that 

oj s y^^ 5 s ...S yjj . 

In view of Eq. (20) , this implies that for any eigenvalue X of A, 

2^.2 ,2^2 
< -X or X < 

meaning X is purely imaginary and |X| > (J^^. Thus, if we denote the eigenvalues 
of A by Xj , j = 1,2,...,2N, then we may arrange them as 

Xj ~ i>^yj 1 ~ 1,2,. .,,N 1 

and ( (41) 

X. =-i»^y~" = X* „ j = N+1,...,2N I 
j ^j-N 2-^ ’ ’ / 

Although Eq. (31) was useful in proving the Fact 3, it is not appropriate for 
further localization of the eigenvalues utilizing the noted properties of 
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dj(y). 

because the Gj^(y) In Eq. (32) contains d^(y) (j < A) . 

To avoid this 

difficulty, we may expand the products dj(y) dj^(y) In G^(P) 

and collect all 

the terms Including Eq. (31) to obtain (see Appendix C) 


N 



g(y) = 1 - L p„ d,(y) 
£=1 

(42) 

where 

N 

_ A 1 ZT ^ 1 



+ 52 ‘^j£ ‘Sdl ^jk 

1 S j < k « N 

(43) 


j.k i- 


and 

2 

A 2 

^2 _ ^2 
^ j 

(44) 

with 

djj^ > 0 If j < £ 

(45a) 


dj^ < 0 If j > £ 

(45b) 

2 

Consider the case when y approaches to Then, d^(y) 

grows without 

bound. 

so we have 



11m g(y) = - 11m P»dp(y) 

)TKJ^±0 



1 + 00 If > 0 

( ± oo If < 0 

(46) 
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Immediately from the properties of g(y) of Eqs. (28), (30) and (46) follows 
Fact 4, as Illustrated in Fig. 1. 

Fact 4 . If > 0 for i = 2,3,...,N, then the roots of g(p) = 0 are separated 
2 

by the a.'s as follows 


< Pi < 02 < ^2 < ••• < ^ Vijj (47) 



It is evident that > 0 because d^^^ < 0 for j = 2,3,...,N. However, 

Pj^ (£ i 2) is not necessarily of a definite sign; the sign of depends on 

2 1 

the relations between the a.'s and the A"''s. Loosely speaking, P» tends to 

J ^ 

A Z i 

be negative if and A are uncoupled or moderately coupled (meaning and 

H 

possibly are nonzero) and if the two adjacent frequencies o^ and are 
very close (meaning d^^^ is very large). In such a situation it is unlikely 
that there exists a root of g(p) = 0 between those frequencies. This is 
shown by the following example. 

If we rewrite P^ by 


P = P - P 


(48) 


18 


where 


and 


E d..K- 

\ j=^+l 


£ 

jJl ‘'j 


lSj<k<£ £<j<k<N 


£-1 


= E d k| - E <^l£‘Sc£^1k 

J lsj<£<k<N J J 


(49a) 


(49b) 


then for A 0 


> 0 and Pj^ * 0 


in view of Eqs. (45). (See Fig. Cl of Appendix C.) 

Assume that is so close to ^ as to satisfy 


A 

M = d 


rfl 


r,r+l 2 2 

''r+l - “r 


» 1 


and d^( » |djj^| for any j ^ k. Then, we may approximate d^^^ ^ by -1 

and, neglecting terms not containing M, we have 

N 


■M 


and 


P'’’« (k^, - E . d. ) M 

r \ r+1 o r+l,j jr ' 


j=r+2 

?;:«(+ E ^j,r+l ‘^jr 
j=l 


M 


i+l 


r+1 


» (e '-T>Ur*l)‘' 


- E d. 


M 


(50a) 

(50b) 


(50c) 


(50d) 


19 



From £qs. (37) and (39) and Eqs. (50), 


r+-l r 


r+1 r 


implying that 


« -P 

r+1 r 


2 2 

Hence, for CJ^ < y < < 


g(y) » 1 - 


yp. 


yP 


r+1 


:1 -I 


^*-^r »^-^r+l 


y-o^ y-<J 


r+1 


and 


:'(y) = g(y) “ 


a- 


r+1 


(y^^,)" (y-a^,)' 


If we denote the root of g'(y) = 0 by y^, then 


^m r r+1 

and the local extremal of g(y) is given by 

20 

1 - a ^ • 

r+1 r 

If P >0, then P < 0 and 
r rri 

g(y ) < 0 (see Fig 2a) . 
m 

If P^ <0, then > 0 and 




g(y ) > 0 (see Fig 2b) 
m 

2 2 

In any case, no root exists between 0^ < y < 


— y 

(a) 


► M 

(b) 
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. In what follows, we will Investigate some special cases In which some 

12 ^1 

identifiable sets of the are positive. Suppose that A ,A ,...,A have 
the same direction and they are orthogonal to the others, which necessarily 

lie in the plane normal to any A^ , (j = 1,2,...,N ). Then, the system 

N ^ 

1 1 

eigenvalues associated with A ,...,A are independent of the others and in 
addition 

Kj = 0, forj,k,il=l 

implying 

0 OT 

= A A >0 for £=1,2 . 


Therefore, Fact 4 always applies to this group of modes. For the remaining 

(N-Nj^) modes = 0 but in general 0 for £,j C{N^+1, N^+2,...,N} . 

Furthermore, if the remaining modes are classified into two other groups, 
N^+1 N^+N2 jj 

i.e., A ... A and A ...A , which groups are also orthogonal, then 

the problem is decomposed into three single axis problems, each of which can 

Q jLT 

be treated independently with P^ = A A >0, 

3.2 Eigenvalue Localization by the Minimax Theorem 

2 

Taking advantage of the noted properties of the matrices and a , 
another localization of the eigenvalues of Eq. (21) will be accomplished by 
utilizing a result based on the minimax method [4] as stated below. 


Theorem; If A,B,C are symmetric matrices with the eigenvalues and y j > 

respectively, which are arranged in non- increasing order, and 
if A,B and C are related by 


C = A + B, 

then for any s = 1,2,...,N 


/ / / 
; /. t 


1 1 : 


' / ; /,/ ' / /• 
a +, B., .,s''y; +/B, 1 

s / TJ/ ■ 's s / I / / /., 


■ / 
/ 

/ 

/ 


21 / 


//■ / / 


/ 


/. ,, 


/ 


/ 



This means that If we add B to A, all of the eigenvalues of A are changed by 
an amount which lies between the smallest and greatest of the eigenvalues of 

B. 

In order to apply this theorem, the characteristic equation of Eq. (21) 

Is rewritten without changing the eigenvalues as follows. 

,.-li 


Since I 4 0 and ja| 0, 

lyUjj - M^a^l = 0 


\ 


Is equivalent to 


0 


2| I -li 


or 


- U„ =0 

1 N 


(51) 


If we define V by 

then we may rewrite Eq. (51) In terms of V as 

IvUjj - = 0 

In which the matrix o ^ Is also symmetric and positive definite. 

If we partition the matrices A and o Into submatrices by writing 


(52) 


(53) 


and 


• ■ [• r-] 


(54) 


(55) 


22 



where 


then 


so 


and 





V 


• 


/V 

• 

a = 

• 

• 

and 

A = 

• 

• 


’0 , 


• 


N 

1 

-a\ 


A^l" 1 A 

1^ 


AA^ = 

t 

- — — 



A A^" 1 A 



-1 

. .T 



m/ = 

U - M 



1 

# 




r" 

1 

“1 


, .I.IT I .1;^ 
1 - A A I -A A 


* .IT L, 

-i A Ws_i - M 




o-\i - a^a1")o-i 


— 0 Au 0- 


-<.-1 A ^? 


A_1 -'A- 

C^N-l ■ ^ 


( 56 ) 


(57) 


Writing 




(58) 


where 


and 


/ 1 A 

’ 3)^ = 


-1,, .l.lT. -1 

a (1 - A A )a 

— ^ I 


0 

/ 


1 A 
= 


^ 1 ^ ■! 
o' CVl ■ AA )a~^ 


-1.1;^ '"-1 
-o^ A A a 


- a AA 


we denote the eigenvalues of O S> ^ and by , vj and O)^ , 

respectively. Note that the matrices and are also symmetric. 


(59) 


( 60 ) 
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The eigenvalues, vj, of SI'*’ are the roots of 


Iv'Ujj - 2>^1 = 0 


l.e. , 


v' - 




IM'T, >1 


or 




|v' - o"^(l - iV‘')cj'^|-|v’TJjj_^ - = 0 


Obviously, one of the eigenvalues of to be called vj^. Is given by 

Vi . 0^1(1 - - A” 


(61) 


The eigenvalues of 9 are the roots of 


ulJjj - ^1 = 0 


(62) 


or 

(0 1 



;>1;^.1T -1 ! 

-a AA 1 

%-l 


- 0 


For the nonzero to the determinant Identity [6] applies to yield 
|( 0 Ujj_^| |(o - = 0 


which reduces to 


2 -1.1^'‘-2;;'.1T -1 

to - A A a AA 


= 0 


( 63 ) 


It turns out that there exist only two nonzero roots of Eq. (62), and these 
are equal In magnitude and opposite In sign. If we denote the positive root 
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of Eq. (63) by o)^, then the eigenvalues of are arranged in non-increasing 


order as 




(N-2) zeros 


where 

™ . -1 (f 

1 - 

As stated previously, the matrices 2)^ and are both symmetric so 

that we may apply the minimax theorem to the matrices of Eq. (58) to find 

- 0)^ S S (6 

From Eqs. (61) and (64), 

1 aIaIT J-.IT 

1-AA .,vAAa AA 

V * ± u) =5 ± 

' 1 

1 - aI^ii t V^i^icoyi)' aaI^ 


On the other hand, if the definition of v (Eq. (52)) is substituted into 
Eq. (65) , we have 

1 - o)i ^ ^1 ^ 1 - 00^ 


i- = -^ 

1 v; 


1 .I.IT . ,/71aT, -^1x2 ^.1T 

1-AA + VAA (a^a ) AA 


" V 


■1 aIaIT y.l^T, ^1.2 ;^.1T 
1-AA -VAA. (Oj^a ) AA 
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We could apply the mlnimax theorem to the other roots of O ^ to 

produce a similar result to Eq. (65). However, the eigenvalues of other 

them are not immediately feasible in view of the equations immediately 

preceding Eq. (61). For those roots of which are expected to give the 

bounds of V (s=2,...,N), we can reformulate the eigenvalue problem of 
8 

Eq. (53) as follows. 

Let be an elementary matrix [6] which interchanges the j-th row 

(column) of an arbitrary square matrix A with the k-th row (column) if A is 
pre-multlplled (post-multiplied) by V^. Specifically, takes the form 


= 


! 


u 


I 


"m-1 


K - ^ 
I I 
I I 
I I 
I I 
I I 
I I 
I I 


“N-k 


j-th 


k-th 


(67) 


f 

j-th 


+ 

k-th 


Note that 




C68a) 


(vJ,2 . (vJ)(v3)’^ - (vJ)’'(vJ) . 


Ivjl^ - 1 


(68b) 


If we define V by 


/ = 

2 3 4 •" £ 


(69) 
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then from Eqs. (68) , 


= (V^ V3 ... vJ"b(V2 V3 ... vj 


= V2 V3 ... ... V3 V2 


= U., 


(70a) 


and similarly 
V 


W / . n 


(70b) 


Premultiplying Eq. (53) by and postmultlp lying by_V?^_ yields^ 

|vvV^ - v^a"^ = 0 

From Eqs . ( 70) , we may write 

|vUj^ - (v^a“ V^) (v^m^M”^) (v^o'M'^) I = 0 I 

i i 

By the noted properties of Vj- composing V , It follows that 


vV^ = 


'k 

-1 

0 






!i i 


and 




1 - A A 


^5. J.T 

"»-r ‘ ‘ 


where 


a 

a = 


'l-l 


i+1 




(71) 
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By these definitions, we may reformulate the eigenvalue problem of Eq, (53) 


- C 21 ^ I = 0 


where 


-1,, .£.£T. -1 I 

s . i 1 ; . . i. 








-1 —1 I 

-(a^) ^ AV\a^) ^ j 


for £ = 2,3, ... ,N. 


Following the procedures taken for we have similar results for 
Vj^(s “ 2,3,... ,N) i 

'■i - * '’j ^ '’1 “i 


where 


, ^ -1,, .£.£T, -1 1 - A A 

vj - (1 - 4 A )Cj J- 
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and 


and 





2 



1 - aV" 

2 


1 - A A -\j h A (a /a^) A A 
for «, = 2,3,...,N . 


( 78 ) 


(79) 


Thus, we have obtained a range in which ]Jjj^(£=2 ,3, * . . ,N) exists, and the range 

I ''H 

is easily evaluated by computing the scalar quantity consisting of A , A and 

a as defined by Eqs. (71) and (72). The ranges given by Eq. (79) may be very 

useful for localization of the y if the quantity in the square root is so 

s 

small that we may write 



(80) 


SL 1 

Strict equality holds if A is orthogonal to any A“^ ( j ^ Z) so that 


. 0 


In general, however, the quantity is not necessarily small because some 
elements of the matrix (o /a^) can be very large. In fact, (o /a^ takes 
the form 
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-9 




' 2, 2 






in which the first (£-1) elements are equal to or greater than unity in view 
of the ordering of Eq. (4a) . 

3.3 Multiplicity of Eigenvalues 

We will first examine the case when 's are all distinct. From Eqs. 


(24a) bnd (42), 


'/ 


N 


N 


N 


f(y) = 77 - y E Pj, 77 

m=l i=l m=l 

s 


(81) 


because |yUj^ - 0 vn (y-<j^) . Although Eq. (42) is derived under the 

m=l 2 

assunqjtion (see Eq. (23)) that |yUjj - O | 0, Eq. (81) is not so restricted 

(the restriction can be removed by the continuity argument), and holds for any 
real value of p for the^distinct Oj case. 


We may rewrite f(y)) as 

N 


f(p) = (y-<^s) 


N 


n - w r Pj n <.v-oi) 


m=l 

m?*s 


i=l ~ m=l 
I’^s m^JljS 

-<jI p. 77 

• * s m 


s s 


(82) 


m=l 

m^s 


which indicates that 


f(o‘) . 0 
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if and only if P =0. We have noted that > 0 if iO, so that 

^2 1 ! Z 

cannot be the eigenvalue, of M^o unless A =0. However, even if A ^ 0, 


can be zero for £=2, ... . ,N, so that a(,(£=2, . . . ,N) can be the eigenvalue; of 


Moreover, if P^ = 0, we may write 


f(y) 


= (y-ap 


-J. ^ ii. i 

fj (y-cf3 - y 1C Pjj^ fj (y-a^) 
m=l “ £=1 ^t=l “i 


m— X 
m^s 


in which a is one of the roots of 
s 


N 

1-y E 

jt=i 

S,^s 


2 

y-aj^ 


1 - y 52 


= 0 


A=i y-<J{ 


if 


1 E 


‘A 


A=1 - a, 

A^s " ^ 


= 0 . 


M^O 


(83) 


This means that if P =0 and if Eq. (83) holds, y = Og is a repeated eigen- 


value of M^cr with multiplicity 2 (at least) . 


L 


Consider next the case when is repeated with /multiplicity /2, /i.e.. 


Og = assuming that all others are distinct. Then, Eq. (31)^ may be 

written as 

N ( A-1 


g(y) = 1 - E d.(y) aV^ - E d.(y)R* 

A=1 ^ i=l J ^ 


lSj<kSA7l ^ 


d.(y)d^(y)Lj^ 


/ 
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/ ' 


N 


,s+l 


g(y) = + {d^(y)}^ - Z {dg(y)}^ + s(m) 

^“S^2 * 


and subsequently 


f(u) - i^(w^3[(d (m)F - E '■s,s+lK<>'>}\(>‘> + '(I') 

hrI * Jl=s+2 


/\ 2 

where = 0, or 


N N j 

f(li) - /7 (M-O^) - E Lf dj(M) + «M) 

m=l £=s+2 ’ J 


(84) 


nR^s ,3+1 


It follows that 


f (c^:) = 0 . 


if and only If 


,3+1 


N 


^ ^3 8+1 = ° 

® £=s+2 ^ 


(85) 


2 

This means that if is repeated with multiplicity 2 and if Eq. (85) 

2 

holds, then a is one of the eigenvalues with multiplicity at least unity, 
s 

2 

If the multiplicity of a is three, i.e., a =0,^=0,^ and others 

S S S rt S^^ 

are all distinct, then as previously 

N 

f(y) = - 77 (M-a^) • + f(y) (86) 

m=l ’ 


where 


0^0 
m s 


f(y) = 0 for y = a^, and hence 


f(a„) = 0 


s+1 

if and only if L . , = 0. 

s,s+l 
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2 2 

If the multiplicity, r , of a is greater than three, then a is one of 

s s s , 

2 I 

the eigenvalues of M O , with multiplicity r —3 (at least) . 

X 8 

The results obtained in Section 3.3 are summarized as follows; 

If are all distinct and of Eq. (43) are all positive for 

2 

£.=1 ,2, , . . ,N, then the system eigenvalues of M a are all distinct as stated 

1 

in Fact 4. , ! 

If CJp's are all distinct but P =0 for some s, then O is the eigenvalue 


of with multiplicity at least unity. 

2 

If a 's are repeated with multiplicity 2, then O is an eigenvalue of 
s s i 

2 

if Eq. (85) holds. Otherwise, it is not so. 

2 

If 0 's are repeated with multiplicity 3, then 0 is an eigenvalue of 
s s 

2 

Mj^0 if Eq. (86) holds. Otherwise, it is not so. 

2 

If are repeated with multiplicity r^ for r^ > 3, then 0^ is an 

2 

eigenvalue of M,0 with multiplicity at least (r -3) . 

X S i . . 
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4. EIGENVECTORS 


Consider the eigenvector problem of Eq. (21) premultiplied by M. 


,-l 


= 0 . (87) 

Since (= - AA ) is symmetric and positive definite and so la^ a'^, the 

results in [1] and [4] apply directly to this problem to yield: 


< 12 N 

Fact 5 . There exist N independent eigenvectors to be called \p ,\{f 

» , 

regardless of multiplicity of the eigenvalues. Furthermore, we may choose a 

set of ijP -s such that and are orthogonal with respect to i.e., 

= 0 for j k . 

'Ik 2 

I whereby and i/; are also orthogonal with respect to a , I.e. 

j 

= 0 for j k 
If we normalize by 

“1, j = 1,2,...,N, (88a) 

^ jthen 

= j = l,2,..,,N (88b) 

These relations are conveniently expressed in matrix form: if we define NXN 

i 

2 

matrices 'F and A by 

¥ = C89a) 


and 






(89b) 


then 


T -1 

'i' M t = U 
1 N 


(90a) 
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and 

- -A^ 


(90b) 


and 

= -'VA^ (90c) 

In what follows, it will be shown that the particular structure of M~^ 
simplifies the eigenvector calculations. 

Assume again, as in Eq. (23), that |v>jU^ “ o | ^0 for Jij . Then, from 
Eq. (87) with Eq. (6e) , 

Vij(Ujj - M'*')ijr^ - = 0 

so that 

(MjUjj - - Pj mV = 0 


or 

(91) 

It is noted that aV is a 3x1 matrix so that if we define by 

= PjAV , (92a) 

then we may write as 

= (PjUjj - o^)'^Ax’ (92b) 

/ 

y 

Subsequently, the eigenproblem of Eq. (87) becomes 
0 = (jjjUjj - a^)i/P - PjAaV 

= MjAx^ - y^M^^ (PjUjj - a^)"^Ax^ 

= PjA{U^ - PjA^ (y^Ujjj - a^)“^A}x^ 
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(93) 


On the other hand, the characteristic equation of reduced form (Eq. (24b)) 
assures that the determinant of the matrix In the braces is zero for , and 
hence there exists a non-trivldl solution for determined by 

{U3 - - 0^)"^A}x^ = 0 

in which only two elements of are to be solved. The normalization con- 
dition of Eq. (88a) with Eq. (92b) substituted yields 

X^V(M^U^ - a^)"^(U^ - Aa'^)(u^Ujj - - 1 


or in view of Eq. (93) , 


. 2 ,- 2 . 


X-jA‘(UjU^ - ^ ^ — S 


= 1 


(94) 


Thus, we can determine X^(j = 1,2,...,N) corresponding to by Eqs. (93) 
and (94) provided that ^ then the ijp's can be calculated 

from Eq. (92b) . 

Next, we will show that even if |pjU^ - 0^1 =0, the Jp^'a can be 

determined in a similar manner. 

2 2 2 
Suppose that and Pj^ 0^ (j ^ s) ; then - 0 | = 0. In this 

case, we may rewrite Eq. (87) (simply by changing the rows) as 


r \ ' ^1 - j 

• • • • 


1 

■€- 

• M ?r 

1 

'.N.IT '. 2 • .N.NT i .N.ST 

- V ^ • • • • \ - "n - V * 1 -“k* * 


,k 



— - — 

L -\aV" -P^A^A«^ 


L 


0 
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With the definitions of Eq. (71) and (72) , we may write 




aS?ST 1 C^.ST 

- (0 ) - a I 4 

I 


aS^ST 


aS.ST 


1 A 

’i'o 




where ij; is the (N-l)xl matrix. Hence 


/ /C8\2 aS^ST, '>k *S.ST ,k „ 


(95a) 


and 


.S^fST -7k ^ .SaST ,k - 
AA ip+AA il) = 0 
8 


(95b) 


Assuming that A^ 0, we have from Eq. (95b), 




, ■S;;‘ST 

,k A A ,7k 

4'_ s~^ 

® A A 


(96a) 


Substituting Eq. (96a) into Eq. (95a) provides 

I 


\»»-A-<’V-a,A^(»3-^)a-1?.0 


A"A 


's 2 1 

in which - (.o') | 0, so that we can follow the procedures established 

above, simply changing the dimension of the matrices to (N-1) . It immediately 


follows that ij) is generated by 


7k 




(96b) 
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irtiere! x*^ is obtained by solving 

' . / I 




with 




X = 0 (97a) 


or 


the normalization condition 


''WT I ' Ir 

i'p 11 o 

I ! ® 


„kT^ST 
X A 


r* 1 “ 


r 1 

1 0 


<^-B- 

0 1 - A^A^^ 



1 


- s J 


= 1 


(97b) 


A A + 


{Wi- (%-i- 


(97c) 


Thus, Vf Cg = andj I ~ 0» is determined by Eqs. (97a) and 

i / 1 ' I i*. ■■ • i' '' V k ■ It 

(97c) , and then i|; is calculated by (96b) and \p by (96a) , so that is 

s 

obtained. Once again, the eigenvector calculation of an (N»l)xl matrix has 
been reduced to simply that of 3x1 matrix. 
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5. TRUNCATION EFFECT 




In this section, the eigenvalues are examined for a truncated system of 
equations represented by 

Y - n = 0 (98a) 


and 


i: -2— — *• 

n + 0ri-AY = O 


where 


and 



(98b) 


(99) 


( 100 ) 


where N is the dimension of the truncated system and 


N < N 


or usually 


N « N 


When such a truncation is exercised, it is required (or sometimes pre- 
sumed) that the nonzero eigenvalues of Eqs. (98) , to be called * 

should be a good approximation of the first 2N smallest modulus eigenvalues 
of the original system, namely A^,A 2 ,...,A^ and Aj^,A 2 ... A^ , as of Eqs. (3) 
or Eq. (8) . 
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(101) 


If we define yj by 

U- ^ 

then the problem is to compare y and y* where y' is a root of 

J *1 J 


and 


|yUg - = 0 


-_T -1 

= (Uj5 - AA ) 


( 102 ) 


(103) 


For convenience of later discussions, we also define "K, F and for the 
part of the system to be deleted by truncation: 


= A 
A = 


*N+1 


LA*^ J 


(104) 


and 


= A 

a = 


*^N+l 






(105) 


(106) 


where 


N = N - N 


(107) 


The notation y! will be used also for the roots of 
J 


|yU - M, 0^1 = 0 
N 


with 


(108) 


j = N + 1, N + 2,...,N. 
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With these definitions , we may write 



In what follows, we will examine the reciprocals of the roots defined by 


. f 1 


as in Eq. (53), which satisfy 


Iv'u- - a 1=0 


vj - l,2,a..,N 


, , =^-l =^1 =^1 , 

|v*U - a M a 1 = 0 

N 


Vj = N + 1, N + 2,...,N. 


Since is partitioned as 


"n I ° 



= _T 
-A A 


I = =iT' 

!U - A A 




- 

-A A 
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the matrix cr ^ is decomposed and rewritten as 


= a) + ^ 


where 


3) = 


_-l i 

a a I 0 

I =i~l=f-l =rf-l 

0 I a a 


(114) 


and 


A 

0 

-1 J ^-1 

-o A fi a 


_T_-1 

-g A A a 

0 

m 


(115) 


It is recognized that the eigenvalues of S> are vj , which are determined by 
Eqs. (112) and (113), because 


, - ' , , _-l_-l_-l 

|vUj^ - S) I = |vu^ - a a 


=-l=^L.-l 
|vU - a M o I 
N ^ 


(116) 


From Eq. (115) , if S =0 then ^ = 0 and 


V. = v! 
J J 


and hence 


4 k.X “• 

In other words, if A =0 for j = 1,2,...,N and k = N + 1,...,N, then the 
eigenvalues of the truncated system together with those of the deleted system 


are identical to the eigenvalues of the original system. 

_ = 1 . 

If A A 0, then the matrix ^ usually causes a discrepancy between V 

and vj . In order to examine this influence, first consider the eigenvalues 

of ^ , to be called o).. The u). are the roots of 

J J 


j 
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0) U„ - ^1=0 
N ' 


or from Eq, (115) 




=-l = _T 

-0 A A a 


_-l_ = =^1 
-a A A a 


uU 


1 


= 0 


For nonzero O), the determinant identity [6] applies to yield 


=-l = _T_-1 _i_-l_ =a? =^1 

(i)U - a A A a (oiUs) a A A o 

N 


= 0 


or 


or 



=“^Q. =T=“^ 

a AA a A A a 


0 


(117) 


(118) 


2 _T_-2_ =0W-2= 

OJ- - (A 0 A) (A 0 A ) 


= 0 


It turns out from Eq. (118) that there exist at most six nonzero eigen- 
values of ^ . Moreover, since ^ is symmetric all the eigenvalues are real, 
so that if we call the roots of Eq. (118) ± u)°, ± a)°» ± ^ 2 

U)° 2: 0, then the are given in non-increasing order by 

OOO^^ o o 

^1 * ^2 * ^3 * y 

(N-6) zeros 


From the fact that the matrix 0 ^ is a sum of two matrices, ® and 

, which are both symmetric, we may apply the Wielandt-Hof fman theorem [4], 


[8] stating that if and vj and u)j are arranged in non-increasing (or non- 

decreasing) order, then 

N N 

H (v, - v')^s 12 (119) 


j=l 


j 


j=l 
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Noting that 


E 0)^ = 2 E (0)°)^ 

j=l J j=l ^ 


= 2 trCS^o A)(^ o ^) , 


and recalling the definitions of Eqs. (52) and (111) , we see that Eq. (119) 


Immediately Implies that 
N 


Hi --i) 


2 _T _-2 -T^ -2= 

s 2 tr(S a h)(E a A ) 


( 120 ) 


Eq. (120) gives a bound of errors which could take place if we approximated 
'|ij by u j . If the relative errors are small, l.e., 

IV'il 1 


where, 


-then 




1 1 


1 1 1 £ ^ ^ , 2. 




Therefore, 


z:(^- 4 y = Euy(^)\oie 

3=1 V^j j=l\»^£/ 


and if we neglect 0(££^) > then 


N 


j 




Each term in the siommation of the left hand side is positive so that for any I 
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or 



_T _-2 _ jr 

tr.(S a 5) (A a A ) 






tr. (S 


l)Ct 


A) 


il = 1,2,... ,N 


( 121 ) 


Equation (121) gives a bound for the relative errors, and if the dimen- 
sionless quantity of the right-hand side of Eq. (121) is sufficiently small, 
judging from some practical point of view, then we may employ as a satis- 
factory approximation of , permitting the truncation to be acceptable. 

This condition, however, is a sufficient one in that even if the quantity in 
question is not small enough this does not necessarily imply that the yj are 
unacceptable. This is because the error limits of Eq. (121) are overestimated 
due to the neglect of the positive terms in the left-hand side of Eq. (121). 

Equation (120) also implies that 
No 

L / 1 IV Tv 

( TT - 1 ^ 2 tr. (S 0 S)(A 0 A) (122) 

j=l \^j ^j / 


The requirement that the quantity of the right hand side of Eq. (122) is 


practically small- is also a sufficient condition in the sense stated above. 

The minimax theorem as applied to the matrices 0.^M^^0 and 

produces the following result. Since the eigenvalues of 0 ^ and are 

Vj and Vj , respectively, and the maximum and minimum eigenvalues of 9 are 
(i)° and we have 


• o . o 

V - CO- ^ V. ^ vl + 00- 

J 1 J J 1 


(123) 
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or 








1 + p j W° ^ 1 - lljO)® 


(124) 


Equation (124) affords an explicit error bound for pj at the cost of solving 
the eigenvalue problem of ^ , which is a 3^3 symmetric matrix. 



1 
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6. DAMPING EFFECT 


This section considers the effect of damping, represented by D “ 2^0 in 
Eq. (14) , on the system eigenvalues and eigenvectors. Rewriting Eq. (14) as 


A ^2N “ ‘ 


- 

1 

0 



> N 

2 

1 

1 

-M,a 

1 0 

1 

1 


“ 

0 

! . 0 



0 

1 -M.D. 



«0 (125) 


we treat the second matrix Including D as a small perturbation to the first 
matrix in the parentheses, whose eigenvalues and eigenvectors have been 
characterized to some extent in the previous sections. We employ the pertur- 
bation method established for dlagonallzable matrices in [7]. Define 


A 

0 1 
1 




„ 2 1 
L'V 1 

0 



Let e be the maximum damping ratio of the appendage vibration mode, i.e.. 


and define pj by 


then 


e = max C 

j 




(126) 


A 

Pj ” g > j ~ 1»2,...,N 


0 < Pj SS 1 . 


By nature of the damping ratio, Cj is a small positive number (typically 
0.001 ^ ^ j ^ 0.01) and so is £, i.e., 

0 < e « 1 . 
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Defining 



we may rewrite the perturbation matrix as 



where 


„ A r 0 I o' 
■ 6 I ‘ -2Mj^pa 


( 127 ) 


(128) 


For e = 0 (unperturbed system) , the previous results indicate that the 
eigenvalues of are given by 







1,2;. ..;n- 


* 

> 



j = N+l 2 N 



to which correspond the eigenvectors 


Y^, 




j = 1,2,. ..,2N, 


From the orthogonality of ijp together with, the normalization condition of 
Eq. (87), the Y^'s are orthogonal with respect to Kj defined by 


A 


1 M-1 

2 “l 


0 • 


— ^ M ^ 
2 1 
2XT ^ 

■ J 


( 129 ) 


48 



In fact 






2 1 I 


' 2X^ ^ 

I j 


.’id .' 

Xj ijr* 


/v^v +(^)i''‘Xv. 


0 

1 


If k¥ i 
if k - J 


(i30) 


In the presence of nonzero e, we write the eigenvector as Y^(e) corres- 

j ; 

ponding to the eigenvalue' (c) and hence the eigenvalue problem becomes 


I ' 


where 


XjCe) ^(e) j Y^Ce) = 0, j = 1,2,...,2N 


orf(e) = |c«/ + e P 


(131) 


C132) 


For a sufficiently small we may assume that 


Xj(e) “ Xj + + a2^ +... 

Y^(e) = Y^ + e +... 


(133) 


(134) 


1 2 2N — 

Since {Y"'', Y Y'^”} form a complete set of eigenvectors, any vector in the 

2N-dlmenslonal space can be expressed by a linear combination of Y^'s, so 


that we may write 


i • 2N 

k=l 


i = 1,2,. . . 


(135) 


Substituting Eq. (135) into Eq. (134) yields 
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+ e 


2n 


2N 


Y^(e) - 


E 


k=l 


Ik 


Y + e 




k=l 


= (1 + ecjj + C^j +...) Y^ 


+ (e +...) Y^ +... 


+ (e +...) Y^ 


If we normalize the coefficient of Y** to unity. 


Y^(e) = Y^ +( S t. ) yV...+ ( S t 

,k=l / Vk»l 

2N / «> 


,2N 


Z L e--t 

Jl=l \k=l 
A>‘j ^ 


kf, 


kN 


or 


2N 


2M 


Y^(e) = Y^ + e Z t.j Y^ + Z 
A=1 A»1 






t2j^ Y +... 


Substituting Eqs. (133) and (137) into Eq. (131) yields 

(a/ + ep) Iy^ + e Z t,„Y^ + Z +...! 




- (Aj+ +...) 


■11 




■21 


Y^+ e Z t,„Y^ + Z t,„Y^ +... 


'U 


■2V 


Equating the coefficient of e in the left hand side to zero produces 


a/ Z t,.Y^ + PY^ = X, Z t,.Y^ + a, Y^ 


(136) 


(137) 


= 0 


(138) 
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In view of the unperturbed relations 

^ ^ (i = l,2,w..,2N) 


we have 


^ E t Y^ = 23 X t ^Y^ 


and Eq. (138) becomes 




(139) 


Determination of 

1T 

Premultiplying both sides of Eq. (139) by Y"* and considering the 
orthonormallty of Eq, (130), we have 

Kj P . 

From Eqs. (128) and (129) , 



- 1 



I “iM ® 


0- ! 0 . 

K- P = 

0 

1 2Xj ^ 

r 


0 1 -2Mj^po 


0 1 0 


8 

• 1 

0 1 — s- pa 

L 1 J 



Hence 



I j I 

t 

1 -1 
1 

t\ . n 



“l = 


U 1 U 

1 “X - 

0 l - pa 

L ' ^ J 


- F_. «. 
xV 

L J ~ J 




(140) 
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and 


a^e “ • e 


(141) 


Therefore, 


XyCe) s'Aj - I 


(142) 


Since is real and positive, the perturbed eigenvalue ^j(e) has a negative 

real part. Moreover, the eigenfrequency |Xj| Is. not perturbed to within the 
first order of e. ; 

DetemlnaClon of t 


Ik 


kT 

Premultiplying both sides of Eq. (139) by Y with k J , we have 
PY^ + (\-^j)tj,jj Y^Ky Y*^ => 0 , 


In which 


and 


YkTR^ Y*" - 1 


Y*^K^ PyJ 



T ■ 

. . 



I 1 

I “ 2 

1 < 


pa 











Subsequently, 




k = 1,2,...,2N; 
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k ^ j 


or 


^ . 


•Ik - \ ^ 


,kT. V 

ijj poijr^ ... 


j 


k =» 1,2,. ..,2N 




Therefore, the coefficient of e In Y-'(£) of Eq. (147) becomes 

X. 


E t 


5 


.and 


•Y^ (e) w y 4 .+ e E ( ’^) x^^ - V ' ^ 

It should be noted that Eq. (143) is valid only when X^^ and X^ are 


for k j so that tj^j^ is sufficiently small. 


/' 
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\ ' 



;(1A3) 

apart 



7. CMONICAL TRANSFOHMATION AND OONTKOLLABILITY EVALUATION 


This section discusses how to construct a transformation which carries 
the system equations (Eqs. C5)) Into a canonical form, and provides Inter- 
pretations of system controllability In terms of the transformed equations. 
We consider the transformation 


I 



(144) 


where T la a (2N+6) x(2N+6) matrix and Z Is a (2N+6)xi matrix, representing the 
system vibration modes. The matrix T Is conveniently partitioned Into the 
submatrices 


T 



(145) 


where ^ 2 ^^ and T 22 are 6x2N, 2Nx6 and 2Nx2N matrices, 

respectively, and the eigenvectors and generalized eigenvectors corresponding 
to X=0 may be assigned to the first six columns, and the eigenvectors cor- 
responding to the nonzero eigenvalues to the remaining 2N columns. Then, 
from Eqs. (15) and (16), we may write 


Til = 


(2) „(2) ^(3) _(3) 


g 


g 


"1 0 0 0 0 O' 

0 0 1 0 0 0 

0 0 0 0 1 0 

0 1 0 0 0 0 

0 0 0 1 0 0 

.0 0 0 0 ff' 1 . 

/ 

/ 


(146) 
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and 


T 


21 


0 . 


In view of the identity (not the transformation) 


'^1 


1 

• 

1 

^2 


Yl 

^3 


Y 2 

• 

= T,, 

• 

^1 

11 

^2 

'*'2 


^3 

• 


• 

3 J 


1^3 J 


(1A7) 


we may well write the variables of Z corresponding to X = 0 as the three 2x1 
partitions 






j = 1,2,3 


(148) 


4 5 

and the two Nxl partitions 2 and Z , so that 


T 

L Z^ J 


(149) 


4 5 

where Z and Z represent the contribution of appendage vibration modes of 
number 2N with nonzero frequencies. 

2 2 

For the nonzero eigenvalues which are determined by |M^a + X = 0, 

we have from Eq. (41) 


X. = i 
1 J 


and 


j — 1,2,.,.,N 
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We know that i|r^ satisfying Eq. (21) is the eigenvector (of x^) corresponding 
to Xj and the zero damping case. In view of Eqs. (17) with D = 0, 

we may write the eigenvectors of the system corresponding to the X^ in terms 
of Xj , and the parameter matrices as 

1 1 T 2 -i 

= - -T A M.a (150a) 

x: 


x^ = \J;^ 


*4 “ 


j = 1,2,... ,N . 


If we define A and A by 
^1 

A = ^2 


and A = 


then Eqs. (150) are expressed in matrix form ((2N+6)x2N) 
r 1 .2 N ' N+1 . 2N~| r .T„ ^2, .,,-2 ' , 


2 

X, — • 

N 

•X, 

1 

1 

N+1 

X, - 

2N 

— — X, 

1 

1 

1 

1 

1 

2 

x« — - 

N 

•”X« 

1 

1 

N+1 
x« — 

2N 

— -X-, 

2 

2 

1 

2 

2 

2 

— — 

N 

1 

1 

N+1 

- 

2N 

— — x^ 

3 

3 

1 

3 

3 

2 

N 

1 

N+1 

2N 

X 4 - - 

"*4 

1 

1 

^4 - 

- 


T 2 —2 

-A^n^a 'HA 

T 2—1 
-A^Mj^a fA ^ 


T 2 *— 2 

-A^M^O 'i'A 

T 2 *-l 

-A^M^a M 


(152) 
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\ 


The upper and lower submatrices of Eq. Cl52) are Identified as ^22 of 

Eq. (145), respectively, l.e.. 


and 


T 2—2 

-A^M^a M 

p 

1 

1 

T 2 *— 2 

-A^Mj^O 'FA 

T 2-1 

-A M^0 VA ^ 

-|- 

1 

T 2 *— 1 

-A^M^O 'FA ^ 



(153) 


(154) 


Thus, we have 
given by Eqs. 
Since It 


established the transformation matrix, T, whose submatrices are 


(146), (147), (153) and (154). 

^^1 = 1 0, and IT22I 0 fey Fact 5, and since T2^^ = 0, 


T 




T Is also nonsingular and 


t“^ = 

1 

-T~ T 

11 12 22 


0 } 
1 

T~^ 

22 


where 


and 


-1 T 
T = T 
11 11 


1 
0 
0 
0 
0 
L 0 


0 

0 

1 

0 

0 

0 


0 

0 

0 

0 

1 

0 


0 

1 

0 

0 

0 

0 


-1^1 r I 

22 “ 2 [- --1 I -*:i--lj 


0 

0 

0 

1 

0 

0 


0- 

0 

0 

0 

0 

1 


(155) 


(156) 


(157) 


. (158) 


Next, we will show that the transformation of Eq. (144) transforms Eq. (5) 
Into a canonical form 
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t"^at Z + 


where 


with 


t“^at 


. J 'O' 0__ I 

' “j To 'i 0 

xroTji'u 

o I A lo 

I 0 lA 


(159) 


[::] 


(160) 


Equation (159) Is derived as follows: Partitioning A Into submatrices as 


A = 


Al_| ^12 _ 


0 i A. 


22 


where 


li 




and 



yields 


T"^AT 


T~^ l_T~^ T~^ 

11 1 11 12 22 


. *^11 I ^2 


T 

11 

1 "121 

! -I 


** * ^ “ “^ “ 


M B 

1 

L ° 1 ^22 J 


1 ° ' ^22 _ 


1 

o 

1 ^22 J 


T"^A 
11 11 


^ 11 ^ 12 "^ 11 ^ 12 ^ 22^2 


T~^A 
22 22 


T I T 
111 12 


r 

O' I T 


22 J 


t“1a t I T~^A T + T~^A T - T~^ T~^A T 

_11 ll_ll ' 11_11 12 _ _11 12_22 _ _11 12_22 22_22 

T 

0 I T ■‘"A T 

^ I 22 22 22 / 
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■10 0 
0 0 0 
0 10 


0 0 0 

0 0 1 

.0 0 0 


0 10 
0 0 0 
0 0 1 

1 0 o' 
0 0 0 
0 10 
0 0 0 
0 0 1 
0 0 0 


'o llo olo 0 
0 olo O'O 0 
D"o|oTlo"o 
0 O'O olo 0 
"d~0[7l OjD"! 
0 o|o O'O 0 


0 0 0 
10 0 
0 0 0 


10 0 
0 10 
0 0 1 


1 0 OjO 0 0 

0 0 l|o 0 0 

0 0 olo 1 0 

1 

0 1 olo 0 0 j 

0 0 0|1 0 0 

0 0 OlO 0 iJ 


1 0 0[0 0 0 
0 0 lio 0 0 
0 0 0|0 1 0 
0 1 OlO 0 0 
0 0 o}l 0 0 
0 0 0|0 0 1 . 


I I 

J » 0 1 0 

- - -f— — 

0 I J lO 

b'l'dyj 


T~^A T 
22^22 22 


= 1 J 1 

^ 1 J '0 J L ’I'A ' WA* J 


+ A + A* 

-A*"^'l'"^M^a^4' + A l-A*"^'l'“^M^a^'i' + A* 


But = -4'A^ and A + A* = 0, so 


-1a T =M^-1 • 

'222^22 Lo1a*J 

r 0 ;u^ 1 -A^M^a^TA"^ 
L 0 lo J -A^M^a^fA"^ 


^11^12 


T 2 *-2 

-A M^O^TA 


T 2 *-l 

-A M^O <FA ^ 


T 2 -1 ' T 2 .*-1 

-A M.a^'l'A ^ i-aVo M 
1 I 1 
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/ 


^12^22 


I A 


T 9 * 

-AM, a j 0 


L-AM,. 
0 


A J J. - 

j 'I'A* 


I 

I 

4- 

I 




T T~^A T 
12 22^22 22 


T 2 .-2 ! T 2 *-2 

-A M,a M I -A^-Mj^a M 


-A^Mj^a^VA"^ i" -a'^Mj^o 


2 . *-l 
^'I'A 


I 


A I 0 

" “I 
0 j A 


Hence 


-A^M^a^¥A"^ I -A'^M^a^'i’A*'^ 
-A^M^a^T [ -A^M^a^’l' 


?^ 1^12 '*■ ^ 12'^22 “ ”^ 12 '^ 22 ^ 22'^22 “ ° 


and Eq. (159) is proven. 

■ I ■ . ■ 

Next, we will express T ^ of Eq. (158) in an explicit form. 


-1 ! _-i -1 


0 

11 rhi 12^22 


«2 

0 i T~^ 


0 

L ( 22 J 

1 

L “4 J 


T ~^ [--- 1 - T~^T T ~^ r ~“ 

11 I B 2 I 11 12 22 I 

^22 iV] 

--if o'], if f"'- i A'V'- If 0 "I 

“ r VJ ' 2 L' 'J 

^ T * 1/2 '• 

where A = A + (J) i?(I) 
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-iFol 1 r-A\aVA"M-A\aW'2 r^-ll 
12^22 [-bJ- 2 [-^*-iJ ’■ V 


1 

2 


T 9 — A.*) 

-A M^a '{'(A + A ) 

T 9 —9 A~9 

A Mj^a fCA + A ) 


'i'"^M,A 
1 c 


But A“^ + A*“^ = 0 and A“^ + A*“^ = 2A~^ 


. . T^2'^22 


2 [■ v] ■ ■ [a\Aa- 2] '' 


and with = -4'A^, this becomes 


Hence 


■'l2’^22 [- i J ' [bI(„2)A-2] ■ [aT«^4^ ] 

/ 

PI.. T-P°]J° I 1 ' 

L-B J V22 Lb.J Lb2-^\^J 


But 


Hence 


T T T * T T * 

®2 ” ^ “l^c = ^ ^ ^ + <t> se^) 

= M2 - a"^m^a 
= M 2 - a'*^AM2 
= CU3 - a'^A)M 2 = U3 
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Thus, we have obtained a canonical fom 


.4 

Z 

.5 

Z 




1 . . -\ 

^ • ‘ 

J 

0 

0 


"0“ 

1 0 1 0 





_0_ 

1 0 1 J 1 

1 1 

• A 1 0 

0 

|0|A 


1*. 



4[4'V. 


u (161) 


We can examine the controllability, of the system described by Eq. (5),-^ 

based upon the canonical transformation of Eq. (161). 

12 3 

It is obvious that the modes Z ,.Z .^and Z are all controllable by the 
input u in such a fashion that Z^ is controllable by u. for j=l,2,3 independent 

l..' V J 

of each other. In other words zero frequency modes are controllable by u 

Independently, where u is "normalized" by the definition of Eq. (2c). (Note 

/i ■- , ' ■ 

that the physical implementation of torquers T does not necessarily apply 

' ' ’ i \ 1 . ‘ ' - . ' ^ ■ . 

independently! to each of three! attitude directions.) 

' " ' ■ ■ " 

For nonzero frequency modes, it will suffice to examine the cohtroll- 

■ 4 !! ■ 5 : - 

ability of |Z ^ because it implies that of Z and vice versa as seen from 

^ I ' 

Eq. (161). From Eq. (161), 


.4 ,1 ' 1 4 1 -1-1 

Z^ = A' Z- + 4 A ¥ M,A u 

111 I ' . i c 

■ 1 i 

1 I 
1 I ' . ' 

With the normalization of Eq. (91), 

■ i ■ 

—1 T 

'1' i = r 

so that Eq. (162) becomes 


(162) 


Z^ = A Z^ + ^ A"^ u 

2 c 


(163) 
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I 


-1 ' ■ ! ■ : 

Since A Is diagonal, we may evaluate the controllability of each mode 


based on the presence or absence of nonzero values In the cofrespo: 


Tiding row 


of the matrix ¥ and thereby establish the system coiitrolltolllfcy cphditiohs, 


after noting the multiplicity of the comprising the notizerp e 


A. 


It should be noted that the controllability arguments stated 

A 5 I 

referred to Z (or Z ) but not to n. Subsequently, justlflcat j^on 


based upon the controllability properties Is also referred to Z^, 


lements of 


[above are 

for truncatloii 

e.g.,jif 


L 

it is als( 


• T 4 

I 1 4' A II =0 then Z^ Is uncontrollable and “truncatable" when 
Sc s 

unobservable. However, this does not necessarily Imply [ that n' i^ truncatable 

/® ' 

or uncontrollable in view of the transformation of Eq. (156 ) . i 

At this point, Eq. (162) is to be compared with Eq. (57) of 

says 


[3] , which 


V = a M, V +• A u 
1 c 


If we transform v into w by 


f w 


then Eq. (164) becomes 


or 


m7^ T w = oV(mT^ 4* w) + a u 
1 11 c 


w = f w +• 4’"^Hj^A^u 


But from Eq. (91) , 


4'"^M^ = 4^ 


■ \ ■■ \ 


(165) 


and from Eq. (93) 


4'"HIj^0^4' = - A^ 


I , I 
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so 


w = -A^ w + u 

c 


(166) 


By comparing Eq. (165) with Eq. (163), we recognize that the mode controll- 
ability conditions are Identical, both being determined by the corresponding 
T 

row of 4' A^.- In addition, the system controllability conditions are identical, 

2 

because the multiplicity of the elements In -A is the same as that of A. 

Although Eq. (164) was derived In [3] by rank calculations without 
utilizing any knowledge about the system eigenvalues or eigenvectors, it 
turns out that the Implications of both approaches are identical as far as 
controllability is concerned. 

T 

In either way, the canonical form of Eq. (162) or (163) with T A^ specified 

will be required for further arguments. It should be noted that uncontroll- 
' 4 

ability of w. (or Z ) does not necessarily imply that of any particular 
' s s 

\ • 

element of v (e.g. v ) because of the transform relation of Eq. (165). 

• • ' ■ S . 

4 

However, in some special cases, the controllability condition for Z 

(or w ) does coincide with that of ti (or v ) , and hence truncation of n is 
S ^ 3 S S 

justified from the controllability point of view. For example, this is the 

o i ' ' 

case when = 0 for j s, permitting 



j=l 


*1 "j 


ij)® w 
^s s 


in which 


V 

s 


is affected only by w^ 


and vice-versa. 
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8. CONCLUSION 


i 
/ 

The eigenvalue and eigenvector problem associated with the original 
system of dimension (2N+6) reduces to that of a symmetric and positive definite 
matrix of dimension N with the zero damping assumption (Eq. (21)). The results 

from the analytical method show that the lowest eigenfrequency of the system 

; 

vibration modes is always equal to or greater than the lowesti of the appendage 
vibration frequencies (Fact 3). In some special cases, including the case 
when the system Is decomposable Into three single axis subsystems, the system 
eigenvalues separate the appendage frequencies at least in ^ weak sense (Eq. 

(47) of Fact 4). The results from the mlnlmax characterization localize the 
eigenvalues as given by Eq. (79) . This procedure requires only simple calcu- 
lations of modal matrices. I 

. ■ f 

The multiplicity of the eigenvalues is dependent on the inertial properties 
and modal parameters In a somewhat complicated manner. It is suggested that 
careful examination should be made of the quantities involved, such as 
k!" and L?'. (Eqs. (83)-(86)). 

J ! 

If the appendage natural frequencies are distinct, then expres- 

i « j 

slons for reduced system eigenvalues are always available as follows: 

(a) If = 0, then = 0^; (b) If ^ 0 but = 0, then = 0^; (c) If 
Z 

A ^ 0 and P^ > 0, then the bounds in Eq. (66) apply and ^ in addition, from 
Eq. (47) j ■: 

“l ^“2 “k > 

i 

0 

(d) If A #0 and P^ < 0, then the bounds in Eq. (66) apply. Finally, if the 
appendage natural frequencies are not distinct then, in addition to the 
bounds In Eq. (65), further restricted results are available in Sec. 3.3. 
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The orthogonal properties of the eigenvectors are presented with 
nomalization conditions employed (Eqs. C90)). The (N><1) eigenvectors are 
generated by a 3x1 matrix (Eqs. (92b) and (93)), that reduces the calculations 
greatly. 

Sufficient condition for acceptability of truncation is given by 
Eqs. (120) -(122) as a result of the Wielandt-Hof fman theorem. This evaluation 
is made by the calculation of the trace of a 3x3 matrix (Eq. (120)). An 
explicit error bound is derived in terms of the eigenvalues of the 3x3 matrix 
(Eq. (124)). 

The effect of modal damping is examined by a perturbation method applied 
to the first order form of the eigenvalue problem. The result assures that 
the system eigenvalues have negative real parts and that it does not affect 
the elgei^frequency to within the first order of maximum modal damping (Eq. 
(42)). The eigenvectors can be significantly changed as shown in Eq. (143). 

Based on the eigenvectors of the reduced system, a matrix is constructed 
that transforms the original system equation into the Jordan canonical form, 
which is useful for the system controllability evaluation. The results are 
compared with those of [3] and physical interpretations are given. 

Although 'the exact eigenvalues and eigenvectors are available only by 
numerical 'calculation, the characterizations established in this paper will 

I 

be useful in thi'at they afford some insight into the eigenvalue localization 

1 

and eigenvector properties. These procedures are recommended especially for 

/ 

preliminary analysis, because the requirements for calculations are not 
burdensome — nothing more than algebraic manipulations of matrices such as 

i 

additions and multiplications. 
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Eigenvalue and eigenvector sensitivity analysis is left for further 
research. This is particularly important for the eigenvector characteriza- 
tion, because it may be greatly affected by a small change in the parameters 
(ill-conditioned) , even though the eigenvalues of the system are always well 
conditioned. It is speculated that a perturbation method similar to that 
employed in the damping effect examination will be useful for sensitivity 
analysis, if we can define an appropriate perturbation matrix. 

Since the quantities which give the bounds for the system eigenvalues or 
truncation error are considered to be the Euclidean norm of 3^1 matrices, 
further evaluation of these quantities may be possible by virtue of norm 
characterization (perhaps by linear programming) . 

( 
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APPENDIX A. SOME PBOPEKTIES OF A A AND M, 


1. All the eigenvalues of A A lie between zero and unity. 


Proof. From Eq. (4d), 


with 




^ O ^ . ^1^ 

I - 1 . „ 


(Al) 


where is the Inertia matrix of the appendage about the total vehicle mass 
center c, and q la a matrix representation of the vector from c to the pri- 
mary body mass center, and Is the primary body mass. 


In Eq. (Al) , I^ and qq are both symmetric and positive definite 

A o T 

and hence I - I Is symmetric and positive definite. Therefore, A A is 

/ 

T 

positive definite and the eigenvalues of A A are positive. In addition, 
- A^A = (I*)"^'^^ I°(I*)"^'^^ 

T 

Is also positive definite, so that the eigenvalues of (U^ - A A) , (l-aj)’s, 
are positive. Thus, 


0 < “j < 1 


(A2) 


2. The matrix is positive definite. 

Proof . Let 6^ be the eigenvalues of Then, are the roots of 

I (6 Ujj - Ml) I = 0 


or 


leu^(U^-M^) -0 
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or 




, A e_]^ 

Let Y “ —g — i 

Iy Uj, - aa"^! - 0 

T T 

meaning Y is the eigenvalue of AA . But the eigenvalues of AA are identical 
T 

to those of A A except for (N-3) zero eigenvalues [4] , so that 


and 


1 ■ f. 


withi 


Yj = Oj (j=l,2,3) 


Yj = 0 Cj=4,5 N) 


0 ^ Yj < 1 Cj=l,2,...,N) . 


It follows that 




J-1.2.3 


and 


3 . = 1 


1 4 ,5 , . . . ,N 


implying that 3^ are all positive, j = 
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APPENDIX B. PROPERTIES OF d (y) = 

Assume that the 0^ are all distinct, i.e., 

2 2 2 2 2 2 

- “i> >■ 0 

for 

i ^ j k ^ 1 

Then, 

d^(y)dj(y) = dj^d^(y) + d^jd^Cy) (Bl) 

and 


where 


with 


d,(y)dj(y)d^^(y) = d^^d^^d^Cy) + ^/^jd^Cy) + 


1 > 

d^COj) = 


1 

1-1 



0 

if j > i 

'ij' 

0 

if J < t 


(B2) 


(B3) 


CB4) 


Proof . Under the assumption, we can expand the product of d^Cy) and d^Cy) 
in the partial fraction as follows. 
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\ 




2 

y-a^ 

B' 

1 


2 2 
Ci^j 

\^l 

2 

, "i 

dj^(y) + 

«2_^2 


y o 


1 


2 

y - KT . 


IL 




dj ( y ) 


= dj ^ d ^( y ) + d ^ jdj ( M ) 


Using (Bl), 


d^(y)dj(M)dj^(y) = {dj^d^(y) + d^^dj(y)} dj^(y) 

“ ‘‘’ ‘^Ij 

- <* J 1 f'Sci'‘i<'') ^ * ‘‘ij<\j'“jW), + ■‘3A(''» 


■*Jl’*lk * ■‘t/jk 


2 2 2 2 


IL 


^ 2^2 2^2 
o .- o . o .- a . 
J 1 k j 


« 2^2 
o, —o 


r ^ 2_„2 

'" k -^1 


k ''i L i j 
2 


2 2 2 2 2 2 
a .- a . aj - a ^ a ^- a ^ 


2 2 


2 2 

a, -a. 

k X 


d ^,, d ., 

ik jk 


2 , 2 2 . 2,2 2 . 

^ i (^ k -^ l ) - 

, 2 2 . , 2 2 . 

( ai - ajUak - KTj ) 

/ 2 2 . 2 

( 0 .- 0 . ) o , 

X j k 

/ 2 2 ., 2 2 . 


Thus , we have (B2) . 
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APPENDIX C. DERIVATION OF g(n) , G^(vi) AND P^^ 

(Eqs. (31), (32) AND (43)) 

Associated with the definition of g(y) (Eq. (27)), we define 8 j[^(y) by 

r 2 1 

gj^(y) - 63 - IE dj(y) (Cl) 

where the quantity In the square brackets Indicates the (a,B) element of the 
3 x 3 matrix with a , 6 = 1,2, 3i Note that 


gj,(y) = g(y) 


and 


(C2) 


,(y) = 1 . 


(C3) 


0 0 

If we define 3><3 matrices Q and R by 


and 


then 


.2 ^ 


R 






g;^(y) = IQ 


(C4) 


(C5), 


(C 6 ) 


and 


Z-1 

Q = q - R 


(C7) 


For the determinant of a sum (or difference) of two 3><3 matrices, the follow- 
ing identity holds. 


/| = |Q^-1| _ |R^| + 2 ^ (-l)J+\q.j^|R^j^| -r.j^lQ.j^|) (C8) 

V j *^“1 

A— 1 i 

where q., , r , are the (j ,k) elements of the matrices Q » R , respectively, 

Jk JK 
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£-1 Z 

and Qjj^ and are the (j ,k) minor of Q and R , respectively. (The super- 
scripts on deleted for simplicity of notation.) 

i UTZ 

Noting that R = d^(y)A A , whose rank is at most unity, we have 

|r^|= 0 and |Rjj^( =0 (C9) 

From Eqs. (C6)-(C9) , 


qp(y) “ Iq^I 


|Q^~^ - R^l 






'jk ‘ 


(CIO) 


so that 


J 

gjl(p) = g;i^i(v) - £ (-1 )^'^‘^aJa^1Qj^ 1 dj^(y) 

j »k=l 

Defining (with the superscript on Q , restored) 

J*^ 

G (y) = E (-1 )^'"^aV|q5;^1 , 

J,k-1 J ^ 

we obtain a recursive formula for gj^CV*) ss 


(CIO’) 


(Cll) 


SjO*’ “ *)t-i - 


CC12) 


Note that G^^(y) does not contain d^(y)(j 2 Z) . From Eq. (C12) together with 
Eqs. (C2) and (C3) , we have 


g(y) = 1 - E Gj(y) dp(y) . 

£=1 ^ 

In what follows, G^^(y) of Eq. (Cll) will be expressed more explicitly. 


(C13) 
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For simplicity of notation, define S.i, by 
£-1 


with 


then from Eq. (C4) , 




jk 

i=l 



(C4), 


„£~1 


Q 



and 




Direct expansion of the determinants of the minors [ • produces 


:,( M ) - wj )' {(l - E22 )(i -E33)-(E23)'} 

|(l-E33)(l-E33)-(E33)i 


+ a,) 


“3^5 12^33 (- 233). E33E 

3a5a;|E33(i-E33).E33E 


) 


(C14) 


(C15) 


(C16) 


12 ) 

The terms of the right-hand side of Eq, (16) are conveniently classified as 

(C17) 


Gj,(y) = gJ + cjcy) + G^(y) 


0 12 

where G^(y) and Gj^(y) are the collections of the zero-th, 1-st and 2-nd 
order terms in (°^ ) » respectively. 
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Obviously, 


/ a *\2 ^ ,. 1.2 ^ ,. 1.2 

( Aj ^) + ( A2 ) + ( Aj ) 


- A 


G ^( y ) = -( aJ )^ 


"22 ■*■ ^ 33 i 


-( A ^)^ 


fE ,, ^ E ,, 


/- a ^2 

- ( A2 ) 


K) 


*31 


X^l 


+ 24 ^ A 3 E23 + 2a5aJ S ; 

From Eq. (Cl) 

gJ(w)« -(Ab^ L (Ab^ d,(p) + Z ( AbV ( y)j 
^ ^ (j“l ^ ^ j“l ^ ^ i 

-(Ab^ j £ (Ab^ d,(y) + (Aj)^d (y) 

^ (j=l ^ j=l ^ 

-(Ab^lz (Ab^ d (y) + E (Ab^d (y) 

( j=l ^ 3=1 

.[4 E aJaJ d (yi ^ ., . . 


+ 2A 


Z-1 

+ 2a:a: E 


b^3 


J =1 

i-1 


A^aJ djCy) 


+ 2A^aJ' E Ab| d (y) 
J 1 j»l J i 3 


2,-1 


-E 


d3(u) + CA^i^)^ + CAjA^)^ 

CA*a|)^ + CAjAJ)^ + CA*A^)^ ■ 


- 

- 

2-1 


(CIS) 


= ' E dj(w) {(^4 ■ 4^1^^ ^4^3 ■ 44^^“^ ^44 “ ^i4^' 

j=l J ( 
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If we define by 


Kj = - aJaJ)^ + (A^A^ - aJaJ)^ + (aJaJ - aJa^)^ (C19) 


then It is also written as 


where 


kJ = CA^ A^KA^ A^)"^ 


P 


4 




4 


0 -A- 


j 


(C20) 


(C21) 


It is recognized that 

kJ.o 


(C22) 


with the equality when A^ = 0, in?>lying that the two row matrices have the 
same direction, i^e,, with some scalar a^, 

A^ = a^A*" 


We also recognize that 


(C23) 


P X, X 

With K. thus defined, we may write G-CW as 

J 


il-1 


jJ(li) = - E d (y)K5' 

j=l J J 


CC24) 
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G^(y) = (aJ)^ 

!l e -(e f! 

^ 22^33 \^ 23 l 

f ‘ oh 


'233S33 -K) 

+ (A 3 >^ 

! £332:33 -(Suri 

+ 2A^A^ 

[£332^3 ■ 2^332^; 

+ 2A2A3 

| 2 ; 3 £ 33 - £332^3 

+ 2A3A, . 

1 ^ 32^2 “ ^ 32 ^] 

By the Lagrange's Identity [8], 


2.^3 -(2:33/ 


£-1 


E (A^)^d 

j=l 


Hii: 


(A^3)^d.(p) 


- Z (A^A^ - A^A^)^d.(y)d/y) 
1 Kj<£-1 J 1 J 


Z A^A^d (y) 

£=1 •* J 


and similarly, 

^13^3 " -^12^2 


( £-1 \ /£-! \ /£-! ^ 

Z AiA^3 dj(M|(Z A^^A^3 dj(y)j Z aJa 3 d^(y)j^£ (A^3)2dj(,)^ 

! 

Z CA^A^ - aJ-A^J CA,A^ - A^A^) d. (p) d, Cji) 

l<l<j<£-l ij-ij jzi j 


We have similar relations for the other terms in Eq. CC25) so that 


(C25) 


2 


82 



G,(y) = 




d.(y)dj(y) |(aJ)^A2A^ - A^A^)^ 


■ /A^\2,.1aJ . /A^N^/a^aJ AdAj\' 

+ (A2> (AgA^ - A,A:^) + (A,) (AiA:i - A^A:{) 


13' 


'1“2 


'2“1' 


- aJaJ) 

«a*(a‘a^ - aJa^) • a^(aJaJ - a‘aJ) 
+2A3<aJa^ - A^aJ) • aJ(a‘aJ - A^A^)! 




z 

l<l<j<il-l 


d^(y)dj(y) j aJ(A2A3 - A3AJ) 


+ A^CA^aJ - aJa^) + A^(aJa^ - A^A^) j 


If we designate the square of the quantity In the braces by L^j , then we may 
write 


G^(y) = Z d CWd (y)Lj 


C.C26) 


where 


,.£7l.jT,2 

= (A A A-^ ) 


Since A^(A^A^^) = - A^CA^A^^) = - A^(a5a^^) , we have 


- A li 

^ij Si “ Sj 


It is obvious that 


CC27) 


CC28) 




CC29) 


a "^i iT il i 

with the equality when A (A A"^ ) = 0, which takes place if the two of A ,A 

and A^ have the same direction or if the third yector is perpendicular to tbe 

Z 

cross-product of two vectors. Hence, if Kj = 0 then = 0 (not vice versa). 
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Substituting Eqs. (CIS), (C24) and (C26) into Eq. (C17) , we have 

G.(y) - aV - E K.d (y) + E lJ, d.(y)d.(y) (C30) 

j-1 J » lSi<J£A-l y ^ J 

with the definitions of in Eq. (C19) and of in Eq. (C27) . 

In what follows an alternative form of g(y.) is derived for the distinct 
eigenvalue case. That is, it is shown that if the distinct, 

then g(y) can be expressed in the form of 


N 

g(y) -1- E Pj<i.(y) (C31) 

A-1 * 


z a 

where is a constant expressed in terms of d^^, and previously 
defined. 

We start with collecting terms containing A^(y) In Eq* (C30) rewritten as 


1 - g(y) - .Gj^(y)dj^Cy) +...+ C 4 _iCv)dj^j^(y) 

+ Gj^(y)dj^(y) +...+ G^(y)djj(y) _ 

Since G^(y),... and G^j^(y) do not contain djj^(y) by definition, we are con- 
cerned only with the terms typified by Gj^(y) ‘^Cy) with k 2 £. 

By Eq. (C17) with I replaced by k, 

G^(y)dj^(y) - G^dfcCy) + Gj(y)dj^(y) + G^Cy)d^(y) 

where 

G° dj^(y) - A*"A‘"'^dj^(y) 
and in view of Eq. (Bl) , 
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k-1 

G^(y)dj^(y) - - E dj^(y)dj(y)Kj 

^ j ' 

k-l, . . k-1 ■ 5. . ■ , 


and'^ln view of Eq. (B2) 


G^(M)dj^(y) = E d^(M)dj(y)d^(y) 


l<l<j<k- 


^ t ViAl"} 


= \(M) E 

l<i<j<k-l; ^ ^ .V.,v. :.... .5 i.,, - ^ 

^ Bi<f<k-il 'J j '•« 


For k=Jl, the coefficient of dj^(P) is given by 

A-1 .... 


.£.ZT V* J IT' 
A A - 2^ d,.K 


d.„d,„L: 


^ in 


k=A 


For k > A, the coefficient of d^ y) is 

k-l ■■ • ' ' , ^ 




A-1 


+ ^ ~ £.+l,...,N 


Summing the coefficients of all terms for k > £, 


85 



£-1 


- A A 


N 

+ E 

k=£+l 


“ E d kJ- + Z d d I 
j=l J l<Kj<£-l 


■ij 


k-1 


£-1 


- ^A ^ ^j£‘^k£^£j ^ S 


Changing the dummy Indices for summation Into l,j with 1 < j Is possible, and 

utilizing the Identities (Eqs. (C23) and (C28)) 

£ £ 1 £ i 

^IJ “ Si “ ^£J ^3 = ^£ 

we may rewrite the double summations In Eq. (C32) as follows: 

N . N N 


and 


:=£+! ^ ^ J=£+l J* ^ j=£+l J 


N 


k-1 


El j£l 


£<j<k<N 


'j£‘Sc£Sj 


E S£^i£^£l 

~ ^ S£*^i£Si 

£<Kj<N J*" 

which corresponds to the region In Fig. G1 Implying that in view of Eq. 
(B4) each term Is positive and so Is the summation. 
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Moreover 


N Z-1 


k=£+l 1=1 ^ 




l<i<)l<k<N 


i£“k£"iil 


= Z 

l<i<£,<j<N ij 

which corresponds to the region in Fig. Cl, implying that this summation 
is negative. Finally 


£-1 


P. = A A - 2 ^ 


d.„K: + 




j=l^ J 1 iTjsil 


N 


- Z d,„K; 


d.„d.„L: 


j=r+fJ^ J £<i“sN 


l<i<il<j<N 


(C33) 


We can now identify which term is positive or negative, so if we decompose 


into 


P = P - P 
Jl I 


CC34) 


where P^ and P^ are the collection of terms which are positive and negative, 
respectively, then 

Z-1 


„+ .A. AT V j i^A . V* j j tA 

* - E 


+ Z d d.jL^. 

A<i<j<N 


(C35a) 
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and 


\ = + E ^ 

j=)l+l J ^ l<i<£<j< 


l<i<£<j<N 


In a more compact but less Informative form, P. may be rewritten as 




in 


i,in 


Fig. Cl Illustrates the Implication of the summation In Eq. (C36) . 


(C35b) 


(C36) 
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1 

2 34 

8-1 8 8+1 

N-1 N 

REGION 

SUMMATION LIMITS 

dj<Og ) ) 


© 

1 < i < j < 8 

+ + 

+ 

© 

8 < i < j < N 

- 

+ 

© 

1 < i < 8 < j < N 

+ 

- 


Figure Cl. Implication of 


1 < i < j < N 
ij £ 
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Chapter 2. Determination of Poles and Zeros of Transfer Functions 
for Flexible Spacecraft Attitude Control 

ABSTRACT , A method is presented for determining the poles and zeros of the 

transfer function describing the attitude dynamics of a flexible spacecraft 
characterized by hybrid coordinate equations. It is shown that the problem 
reduces to that of finding the eigenvalues of matrices which are constructed 
by simple manipulations of the inertia and modal parameter matrices. Particu- 
lar emphasis is put on the determination of the zeros, which depend also oh 
the sensor and/or actuator location. The established procedure will be useful 
for numerical determination on the digital computer. 
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i' 

i' 


INTRODUCTION 


The system discussed in Chapter 1 is a multi-variable linear constant 
coefficient system with three inputs and three outputs and a large number 
of states. Such a system is sometimes conveniently described by a "transfer 
function matrix," in which the determination of the poles and zeros is of 
primary importance in designing a conventional compensator by the classical 
approach or an optimal controller by modern techniques. 

The determination of the poles of the system is relatively simple from 
a computational point of view, because it is just a calculation of the system 
eigenvalues, for which computer algorithms and/or programs are well established 
(see also Chapter 1 for eigenvalue characterization) . 

The calculation of the zeros, however, is not an easy task even numeri- 
cally, and even the definition of the "zeros" is not unique. We will first 
survey some of the recent papers dealing with this problem, and thereby 
demonstrate the difficulty of the zero determination problem. 

For the purpose of this survey, we start with a more general description 

* 

of the system (after Laplace transformation of Eq. (5) in Chapter 1 with 
X(0) =0)" 


sX(s) = AX(s) + B uCs) 

with the observation equation (a) 

y(s) = C X(s) 


whose transfer function matrix (sometimes called "frequency response" [1]) is 
given by [1] , [2] 


The Laplace transform of a variable, say X(t) , is written as X(s) , and when 
the distinction between X(t) and X(s) is clear from the context the arguments 
may be omitted. 
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■ G(s) = C CsU„ - A)“H 
so that N 

y(s) = G(s) uCs) 


(b) 


Definition 1; The zeros of the transfer function are those of the individual 
entries in G(s) which give the relationship between the r-th input u^ and the 
j-th output Yj . 

This form of G(s) (sometimes called the first form) arid the definition 
that follows seem to be standard. Rosenbrock, [3] , [4] however, insists that 
the following form (called the second form) is more general; 

Consider the system 


J(s) X(s) = K(s) u(s) 
y(s) = L(s) X(s) + M(s) u(s) 


(c) 


where in particular we may have, as in Eq. (a). 


J(s) = 8 - A, K(s) = B, L(s) = C, and M(s) 

and define the system matrix in polynomial form by 


P(s) = 


J(s)_ ' K(s) 
■-L(8) r M(s) 


which has Smith normal form [4] (uniquely determined) 

0 


N(s) = 


diag(Ej (s)) 


Definition 2 : The zeros of the system are the zeros of the polynomial ^jCs) 

taken all together. 

The formulation by Rosenbrock (Eq. (c)) is certainly more general than 
Eq. (a) since the former permits the case when |t(s) | =0 while the latter 
does not. Usually, the zeros are not identical for the two definitions,* 
except for the single-input, single-output case. 
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In the problem under consideration, |t(s) | ^ 0 so that it is not neces- 
sary to take the form of Eq. (c) . 

In the case where the number of Inputs equals the number of outputs, 
as in the present problem, we may consider the determinant of G(s) , since 
it is square. In this case, the zeros of this determinant are identical to 
the zeros of the definition 2, because the elementary operations tised in 
forming the Smith normal form do not produce excess zeros and the determinant 
remains unchanged except for a scalar multiplier. 

Brockett [2] has presented a way of zero determination by deriving the 

"inverse equation" whose poles Include the zeros of the original system; 

hence the problem reduces to the eigenvalue problem. However, his results 

are limited to the single— input , single-output case, and, as pointed out by 

Davison [5], the formulation of the matrix representing the inverse system 

is computationally difficult, because the prescription for forming the matrix 

relies on logic statements such as a = 0 if d = 0 and a=lifd?*0, 

and because of the finite word length of the digital computer, a subjective 

•*8 “20 

judgment must always be made to determine if d = 10 or d = 10 may be 
taken to equal zero. Also the system must be controllable and observable. 

In addition this approach supplies excess roots at the origin which must be 
sorted out by judgment. 

Davison [5] has proposed a computational method and computer algorithm 
for zero calculation, which also reduces it to eigenvalue calculation of a 
matrix formed by Cramer's method. Although this method also supplies excess 
roots, it is always possible by repeating some steps of the calculation to 
determine which roots should be rejected. The advantage of this method lies 
in the fact that this approach can reduce the inaccuracy due to word length 
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limitation and^ round-off error of the computer. On the other hand, this 
method has a disadvantage for large-scale systems because it requires the 
solution of an eigenproblem of higher dimension than truly necessary, and 
hence the number of excess zeros can be very large. In addition repetition 
of some calculation steps is not desirable. Kropholler and Neale [6] have 
discussed a similar method. 

The method proposed by Guldorzi and Terragni [7] operates only on a 
subsystem of minimal dimension completely describing the dynamical behavior 
between the input and output of Interest. The zeros are finally determined 
as the roots of a polynomial of the minimal dimension (no excess zeros are 

..produced). Although some of the difficulties in Davison's method are over- 

1 

come, f the Guldorzi-Terragnl approach requires the ntimber of Independent 
vectors, which may be dependent on subjective judgment, and the calculation 
of the zeros of a polynomial can be less accurate than that of eigenvalue 
calculation for some problems. 

Yokoyama [8] has discussed a method of obtaining transfer functions by 

transformation to a phase variable canonical form. He has overcome the 

i 

difficulty in the Inversion of CsU - A) in Eq. (b) , but the calculation of 
the carionical form is another burdensome task, so that his method does not 
seem to be appropriate for the present problem. 

In this paper, we will show that the zero determination problem reduces 
to eigenvalue calculations of a matrix. The construction of the matrix is 
discussed in detail for the cases when the sensors and/or the actuators are 
mounted on the primary body. The method presented here does not require any 
repetition of calculation steps or any subjective judgment, unlike most of 
the alternative methods noted. 
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We will consider primarily the zero determination problem in the sense 
of Definition 1. We will also discuss the zeros of the determinant of the 
G(s) which is related to Definition 2. 


DERIVATION OF TRANSFER FUNCTION 


As in Eqs. (1) of Chapter 1, we start- with describing the system in 
terms of hybrid coordinates as 


T.. 

I 0 - 6 fi = T 

fi + Dfi + O^n - 60 = T 


Cla) 

Clb) 


with the observation equation [9] 


where 


and 


y = 0 + 

’ o 

A T 

c c 

A T 


Clc) 

C2a) 

C2b) 


establish the actuator and sensor locations, repsectively, and all other 

matrices are defined as previously. 

For convenience of later discussion, we partition the matrices 6, 

and SB' as 
o 

6 = [a^ a^ a^] C3a) 

i?' = Ib^ b^ b^] (3b) 

c 

and 

SB' = [c^ c^] (3c) 

o 

/ 

where a^ , b^ and c^ (j=l,2,3) are Nxl matrices. 
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Laplace transformation of Eqs. (1) with the zero Initial state (i.e. 
0=0, 6=0, n=0 and n = 0) yields 


2 * 2 T 

s I 0 - s^6 n = T 

(s\, + sD + ahn - 8^60 = 

N C 


(4a) 

(4b) 


y = 0 + i?'^n 

O 


(4c) 


Define the N><N matrix Q(s) by 


Q(s) = s^(s^Ujj + sD + a^) ^ , 


(5) 


then from Eq. (4b) 


n = Q(s)60 + -j Q(s) <?• T 


(6a) 


and from Eq. (4a) with Eq. (6a) substituted 

s^l*0 - 8^6 ^Q(s) 60 = T + 6^ Q(s) .2^ T 


or 


or 


8 ^( 1 * - 6 ”’^Q( 8 ) 6)0 = j «3 + fi’^Q(s) |t 


0 = -| |l* - 6 ’^Q(s) 6 j ^ I U 3 + 6 ^Q(s) j T 


(6b) 


From Eqs. (6a) and (6b) , 


n =- 


Q(s )6 jl* - 6 '^Q( 8 ) 6 j“^ j U 3 + 6 ^Q(s) 


+ Q(s) j?' 

c 


(6c) 


Substituting Eqs. (6b) and (6c) into Eq. (4c) yields 


1 

y = — 

s 


U3 + i?^'^Q(s)6j |l* - 6^Q(s)6j“^ 


\ T 

• ju^ + 6 Q(s)i2’^| 


I T 

' ! + i2^'^Q(s) i?' 
o c 


(7) 
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Thus, we may define the transfer function G(s) from the control Input T to 
the sensor output, excluding the poles at the origin, by 


^ I 


n 


G(s) = {U- + iJ?''^Q(s)6: {I” - 6"Q(s)6 


I'* l"3 


U, + 6'-Q(s) ^ 


+ i?'’^Q(s) i?’ 
o c 


(8a) 


and denote the a-8 element 

of G(s) 

by 


Gii(s) 



G(s) = 

G2i(s) 

0„(s) 

GjjCs) 


G3i(s) 




l.e. , 


(8b) 


Hence, the purpose Is to establish a procedure to detemlne the poles and 


zeros of G 


06 


(a) . 


We will examine the following cases: 


(1) ^' = i?' = 0 

c o 

meaning that the actuators and the sensors are attached to the 
primary body and 


(2) i?’ = 0 and i?* 0 

c o 

meaning that the actuators are on the primary body and the sensors 
on the sub-body characterized by and 

(3) 0 and i? = 0 

c o 

meaning that the actuators are on the sub-body characterized by 
and the sensors on the primary body. 


The reason for omitting the general case when ^ 0 and S' i‘ 0 may 

c o 

may be Justified by the possibility that we may consider either of these sub- 
bodies as the primary body. 
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INVERSION OF (I* - 6 ’^Q(s) 6 ) 


Let F(s) • fF^Ca)] be the adjoint matrix of (I* - 6 ^Q(a) 6 ). Then 


(I* - 6 ’^Q(a) 6 )"^ 


lisi 


I - 6^Q(s)6| 


( 9 ) 


* T 

since, with the definitions of Eqs. (3), I - S Q(s )6 takes the form 


* T 

I - 6 Q(s )6 


^ 1T„ I ' IT- 2 - IT- 3 

I^, - a Qa Ino - a Qa Ino •* a Qa 

2T 1 2T- 2 2T- 3 

I,, - a Qa I,, - a Qa I^^ - a Qa 


'■12 “ -22 

3T„a .3T^_2 


3T^-3 


I, , - a"‘Qa^ 1^3 - a"‘Qa I^, - a''‘Qa 


L"13 


33 


Then the ot -3 element, F^(s) , of FCs) Is expressed by 


F^Cs) - (-l)°^^|l“^ - 6% 5"| a, 3 - 1,2,3 


ao) 


where 


.11 ^ 


.12 A 


- 


- 


r *1 

I I 

^22 23 

22 ^ 

, = 

I I 

11 13 

-33 ^ 
1 1, ® 

^11 ^12 

.^23 ^33- 


.^13 ^3- 


.^12 ^22 

•• • 



, , 

z « 

^12 ^13 

t13 ^ 

» “ 

^12 ^13 

-23 A 
. 1 = 

r t 

^11 ^13 

.^23 ^33- 


I22 I23. 


I I 

_ 12 23 _ 


Cll) 


and 


^21 ^12T -31 U3T ^32 _23T 

I =i , I =r. ,t =1 


and 


J.1 A 2 3, J.2 A 1 3, ^ r X ,2-, 

0 = [a a ] , 0 = {a a J , o = Ja a J 


C12) 


Therefore, we have obtained the Inverse 


Cl* - 


(«1) “^3 ll“^ - ^^QCs)6^| 
1 1* - 6Vs)5| 


C13) 
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It is noted that since I 
symmetric, implying 


6^Q(s)6 


is symmetric. 


the inverse is also 




This is confirmed for example by considering the case of 05=1, 3=2 as follows 


F,2(s) 


12 

- 

Il 3 - 

23 

- 

I33 - 

12 

- 6 ^*^06^ 

I23 - 

13 


I33-6V 


- F 


21 


DETERMINATION OF POLES AND ZEROS 


(1) In the case of = 0, from Eq. (8a) , 


G(s) = (I* - 6^Q(s)6)"^ 


(14a) 


so that, with Eq. (13), the transfer function, G o(s), from the 3-th input to 

ap 


the a-th output becomes 




|I°^ - 6^^0(s)i°l 
I* - 6'^Q(s)6| 


(14b) 


In what follows, it will be shown that the determinants in Eq. (14b) 
reduce to characteristic equations of (2Nx2N) matrices. This reduction comes 
from the observation that the poles of the transfer function must be the 
eigenvalues of the matrix A in Eq. (1) or its equivalent form. 
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Consider the eigenvalues of the matrix of reduced form (Eq. (19a) of 
Chapter 1) 


a/ 


I U 
I N 


- -| 
-Mj^a I 


(15) 


where 


* ~1 T 

= (Ujj - 6(1 ) 6^) 


(16) 


The eigenvalues, s, of are the roots of 




(17) 


But 




sU 


N . 

sUjj+ 


|sUjj|-|sU^ + M^D + M^cy2(sUjj)‘\| 


= IMj^I • ls^M~^ + sD + 

= |R|^|-|8^(Ujj - 6 (I*) ■^6’^) + sD + o^\ 

= 1m^1*|s^I^ + sD + - 8^ 6(L*)“^6’-| 

= + sD + a^|*|Ujj - Q(s)6(I*)“^6’^| 

= + sD + a^l • |l*r^; |l* - 6'^Q(8)6) 


Hence 


1 1* - 6Vs)6| = |M^r^|s^Ujj + sD + 1 1*11 sU^jj - ^1 ( 18 ) 
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similarly, If we define the 2Nx2N matrix by 


■ 06 ^ 




■■ 


I N 


o8 

-M^D 


where 


oB A 
= 




-1 


then the numerator of Eq. Cl4b) Is written 

|i“S-A(«)5“i 

- |rf*®r^|a^U^ + sD + o^f''|l“®| |su 

Substituting Eqs. (18) and (21) Into Eq. (14), we obtain 


2N 


- ji 


oBi 


G^(a) = (-1) 


ori-B 


^ I I s\+sI>+<J^ I • 1 I * I sU,^ - ^ I 


2N 


M^r^|s^Ujj+sI>fa^| • |l’‘| • |sU2jj - 


Noting that 


|M^r^|l*| = ll"! |Ujj - fid") ■•'6 




ii*i-iu. - (i*)”y6| 


I* - 6^6 1 


and similarly 


iM“«r^ii*i = ii“» - 6®^6“i 


we have 




|l“S . |a„^ . ^ 


aBi 


I - 6 6| IsUjj - ^ 


(19) 


(20) 


( 21 ) 


( 22 ) 


(23) 
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Thtis, the poles and the zeros are the eigenvalues of the matrices 
(Eq, (15)) and (Eq. (19)), respectively. The matrix Is con- 

structed by the definitions of Eqs. (11), (12) and (19), but since It contains 
the Inversion of (1°^) and (U„ - J (1°*^) ^6^) Eq. (23) Is restricted to the 
case when 

|l°^| 4 0 (24) 

and 

|Ujj - 6“(l“^)"^ 5^1 y 0 (25) 

However, Eqs. (24) and (25) always hP.ld for the case of a Implying that 

the poles and zeros are always calculated by Eq. (23) If we consider the 
transfer functions from the a-th Input to the a-th output (of=l,2,3) which are 
usually of primary Importance. This comes from the symmetry and positive 
definiteness of the matrices I and - 6(1 ) 6, which guarantee that every 

principal minor of these matrices is also symmetric and positive definite (by 
Sylvester's theorem) with the recognition that is the o3 minor of I and 
|Ujj - 6°^(I*^^)”^ 6 I is related to the aB minor of - (I*) ^6^6. 


(2) In the case of = 0 and 4 0, from Eq. (8a) 
G(s) = {U_ + i?'^Q(s)6} {I* - 6 ^Q(s)6}”^ 


(26) 


By partitioning of Eqs. (3a) and (3c), 

+ J’^^Q(s)6 



r 1T1 




C2T 

‘^3T 
^ c J 

Q(s) [a^ a^ 

a'l 

' 1 + c^'^Qa^ 

IT- 2 
c Qa 

1 

2T- 1 
c Qa 

1 

2T- 2 
+ c Qa 

2T 3 
c Qa 

3T- 1 
L c Qa 


c Qa 1 

, 3T- 3 

+ c Qa _ 


(27) 
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Carrying out the matrix multiplication In Eq. (26) with Eq. (9) substituted 
yields 






E + c°^Qa^)F^.(s) 


I - 6 ^Q 6 | Y=l 


OfY 


W 


(28) 


where 5 denotes the Kronecker delta defined by 

erf ' 


OCY 


1 if a = Y 

0 If a Y 


We record G -(s) for a,B = 1,2,3 as follows, deleting the factors 

ap 

1 1 * - dVi and (- 1 )“^^ 


°ir 


° 12 * 




Si* 


^ 22 * 


S 3 ’ 


Sr 


S2’ 


S 3 ’ 


IT 1 IT 2 IT 3 

(1 + c^Sa^) + c^Sa + c^Sa'' F^^ 


IT 1 IT 2 IT 3 

(1 + c-^^Qa-^) F,, + c-'^Qa'^ F 


32 


IT 1 IT 2 IT 3 

(1 + c^Sa^) F,, + C^^Qa^ F23 + h^Sa-" F33 


12 

7 

13 


+ (1 + C^Sa^) F^^ + c^'^Qa^ F 

c^’^Qa^ F^2 + + c^^Qa^) F 22 + c^Sa^ F 

C^Sa^ F ^2 + (1 + c^Sa^) F 23 + c^^Qa^ F 

c^'^^Qa^ F^^ + c^'^Qa^ F,^ + (1 + c^’^Qa^) F 


c^Sa^ F ^2 c^Sa^ F 22 + (1 + c^SaS F 
+ c^Sa^ F 23 + (1 + c^Sa^) F 


31 

32 

33 

31 

32 

33 


Note: F o = Fq 

oB Ba 


It is easy to verify that Eq, (28) may be rewritten 


i^oB xT„?aBi 

G (s) = (^1)°^^ ^ f 1 . (29) 

^ ll - 6^61 

aqA th. * 

where 1 ^ is the 3x3 matrix formed by replacing the B column of I with the 


3 x 1 matrix 
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and 6 la the Nx3 matrix formed by replacing the 3-th column of 6 with 
(-c^ . We record all the and 6°^ in the following. 

M ^ ^ r*‘ 

^ ^ ^12 ^13 ^,2 ^ ^11 

I = 0 1^2 l23 . I = 1^2 0 1^3 . 1 1,2 

^23 ^33 " ^ 


.V., ^11 ^12 

l"*^"* =1 I 

12 ^22 


0 ^2 ^13 

1 I22 ^23 

0 ^23 ^33 


0 ^12 ^13 

° ^22 ^23 

^ ^23 ^33 


^13 ^23 ° 


/^2, r^ll ^12 ° 
l'-ll,2 I 22 1 

■^13 ^23 ° 


"11 ^12 ° 
^12 ^22 ° 
^13 ^23 ^ 


,12 3 

[-C a a J 


(31a) 


r 1 .,1 3, 

[a -c a J 


(31b) 


r 1 2 1, 

[a a -c J 


(31c) 


r *2 2 3, 

[-C a a J 


(31d) 


r 1 „2 3, 

[a -c a J 


(31e) 


r 1 2 2 , 

[a a -c J 


(31f) 


r 3 2 3, 

[-C a a J 


(31g) 


r 1 J 3, 
[a -c a J 


(31h) 


r 1 2 „3, 

[a a -c J 


(311) 


Eq. (29) is confirmed for a specific combination of indices below; 


If a = 3 = 1, then 
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G^j^(s) *1*1 = (-1)' 


^ ^12 ^13 

° ^22 ^23 

.0 


- 



2T 

“ 

a 


3t 
_ a 


12 3 

Q[-c ,a .a'^] 



IT 1 
1 + a Qc^ 

T IT- 2 

I,, - a Qa 

T 2T„ 2 

^22 “ ® Qa 

' 13 - 

^23 “ 

IT- : 
a Qa 

= 

2T„ 1 
a Qc 

2T- : 
a Qa 


3T^ 1 
a Qc 

T 3T- 2 

I 23 - a Qa 

^33 ~ 

3T- : 
a Qa 


= (1 + a^'^Qc^) + a^'^Qc^ F^2 + a^'^Qc^ F^^ 

= Cl + C^'^Qa^'^) + c^^Qa^ F^^^ + C^^Qa^ F^^ 

= . C-Q aV,.^, 


which is identical to Eq. (28) with a = 3 = 1. 

Recognizing the similarity of Eq. (29) to Eq. (14) , we may rewrite (as 
in Eqs. (21) and (23)) 

|£”«-6^q6“8| 

- |M“®r^|s^U^ + sD + (J^r^ -|l-°*|-|sU2„ - 


so that 


G^(s) = (-1) 


ort-3 


I £“8 - |.U2„ - 


* T 
I - 6 6 


*“2N - 


(32) 


where 




(Uj^ - 6°^(l“^)"^6^) ^ 


(33) 


and 


II > 

1 

1 

1 0 
1 . 

1 

j 

1 

1 

1 

1 1 

L 1 

-m“^d J 


(34) 
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Thus, we have reduced the zero determination problem to the eigenvalue 


,r(K 

problem of the matrix defined by Eq. (34). „ As previously, the construe- 

tion of the matrices and 6^ is straightforward, but since contains 

the inversions of and (Ujj. - 6*^(1^) ^6^), Eq. (32) is restricted to the 


case when 


0 


(35) 


and 




(36) 


In view of Eqs. (30), 

U“*l = |I“®I 


so that Eq. (35) alwaysholds for a' = B. However, Eq. (36) is not neces- 
sarily guaranteed even for a = 0, as opposed to the previous case (Eq. (25)) 
because 6^ contains the sensor location parameter .(-c^ as seen from Eq, (31) . 

It is worth noting that |U^ - 6°®(I°^) ^6^ j in Eq. (36) gives the 

2N ^Ct8 

coefficient of the highest order tern in s, namely s- , and |l^| gives the 
constant term in the numerator polynomial. Therefore, if Eq. (36) is violated, 
then the number of the zeros is less than 2N. In addition if 


|Uj, - 6°^(l“^)"^6''‘^|- |l“^| < 0 


then there exists at least one zero in the right half plane. 

(3) In the case ^ 0 and = 0, in parallel with the preceding case, 
Eq. (8a) becomes 

G(s) = {I* - 6'^Q(s)6r^ {U^ .+ 6 '^Q(s) (37) 
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If we transpose G(s) , then 


G^(s) = {Uj +^^'”^Q(s)6}{l* - 6’’^Q(s)6} (38) 

T 

(Note that Q (s) = Q(s) since it is diagonal.) 

We recognize that Eq. (38) is Identical to Eq. (.26) if is replaced 
by Therefore, we may treat this case simply by changing c with b and 

noting that G^(s) of Eq. (37) is given by Gg^(s) of Eq. (26). All the 
results in case (2) immediately apply to this case. 

ZEROS OF DETERMINANT OF G(s) 


The zeros of |G(s)| are related to those of definition 2 as mentioned 
previously, and they are easily determined for the cases (1), (2) and (3). 
If - Ot then from Eqs. (14a) and (18) with Eq. (16) 


|G(s)| = 


I* - 6’^Q(s)6| 


|s^Uj^ + sD + a^l 


lUj^-6(I*)'V| |l*|lsU2jj- ^1 


(39) 


so that the zeros of | G(s) | satisfy 



+ sD + a 


= 0 


Therefore, 


Sj . * 1 v/TT 






If 


= 0 and St?' 0, then from Eq. 
c o . T ^ 


G(s) 


lU^ + i?^^Q(s)6| 
I* - 6'^Q(s)6| 


(26) 


(40) 


(41) 
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In parallel with the derivation of Eq. (18) 


IU3 +i?^'^Q(s)6| 


U, + 6 I IsV + sD + 0^1 |sU,„ - 


“N 


2N 


(42) 


where 




o A 


0 I U, 


J . N _ 

L -M°a^ ' -M°D - 


(43) 


with 


^ (U - 6 i?’^) 

M O 


-1 


(44) 


Substituting Eqs. (18) and (42) Into Eq. (41) yields 


G(s) 


IU 3 + 5| 

|I* -6^61 


l«^2N - 

®^2N “ ^ 


(45) 


Thus, the zeros are determined by 


IsUjj = 0 


(46) 


which are the eigenvalues of . 

If = 0 and 0, then from Eq. (37) 


|U +6Vs)iZ’:l 

|G(s)|= f- 

|l - 6Vs)6| 


(47) 


which reduces (as in the previous case) to 


G(s) 


U, + 6 ^l\ 


* T 1 
I - 6 6 j 


|sU2^ - ^ 


(48) 
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where 




c A 


0 

I ^ 



_ -N_ . 



1 -m'^d 

_ 


(A9) 


with 


= (U„ - ^’6^) 
M C 


-1 


(50) 


Thus, the zeros of |g(s)| are the eigenvalues of determined by 


latJjN -‘^1=0 


(51) 


The zeros determined by Eq» (40) , Eq. (46) or Eq. (51) do not correspond 
to the zeros of any particular transfer function, » which are determined 

by Eq. (23) or Eq. (32). 

It is pointed out by Brockett [1] that the zeros of |g(s)| play a funda- 
mental role in least-square optimization theory and they are important in 
determining if a "plant" can be decoupled by state variable feedback. 

CONCLUSION ’ 

First we discussed briefly the two different definitions of the zeros 
of the transfer function matrix for linear time-invariant systems, to which 
the attitude control system of flexible spacecraft under consideration belongs. 

The transfer function matrix is obtained for the model of three inputs 
and three outputs (Eq.(8a)). Based on this representation, it is shown th^t 
the zero determination problem reduces to the eigenvalue problem of the 
matrices defined by Eqs. (19) or (34) for the cases when the sensors and/or 
the actuators are on the primary body. These matrices are formed by inter- 
changing the columns or rows of the system matrix (Eq. (15)) with those 
of the sensor or actuator location matrix. They are always feasible for the 


111 



transfer function from the j-th input to the j-th output (j=l,2,3) in the 
primary body instrumentation case. However, for the other cases, the pro- 
cedure is restricted by the assumption of Eqs. C35) and (36). The results 
above stated are related to the zeros of Definition 1, and they will be use- 
ful for designing a controller by means of classical techniques, such as the 
root locus method, 

A method for determining the zeros of the determinant of the transfer 
function matrix is also presented (Eqs. (40), (46) and (51)), and the zeros 
thus determined are related to Definition 2. They are useful for designing 
an optimal controller by modem control theory. 

. . ^ ^Mthough no algorithm has been shown, the procedures may be implemented 
on the, digital Computer, without any repetition of calculation steps or 
any subjective judgment. 
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APPENDIX A. EXAMPLE OF TRANSFER FUNCTION DETERMINATION 

Consider the single axis rotational motion of a spacecraft with a single 
appendage mode as shown In Fig. Al. 

I 



Fig. Al. Model of Example 

The mode shape Is assumed to be given by a function of t.e., (pCO • 
From Eqs. (la), (lb) and (Ic) , with scalar variables and parameters, 


where 


10 - 6 n = T 

.. ' 2 •• 

n + 2?ari + a n - <50 = 0 

y = 0 + 


L 

--f / 


/ 


<J>(C)(d + C)dC 


(Al) 

(A2) 

(A3) 

(A4) 
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T 

and .we Introduce the symbol c to represent , as given by 






(A5) 


Then, the transfer function GCs) from T to y (with the poles at the origin) 
is given by 


G(s) 


1 - 6c| IsU^ - 


111 


I - 6 I \e\J^ - ^ 


_1 

2 


(A6) 


where 





(A7) 


(A8) 


with 


and 


= (1 - 6 ^/ 1 )"^ 


(1 + 6c)~^ 


CA9) 


(AlO) 


Since 


and 


then 


|l - 6^1 ‘laU- \ = l{(l-6^/I)s^ + Ds + 0^} 


|l + 6b|’|sU, - = (1 + 6c)s^ + Ds + , 

^ 1 , 


X i (1 + 6c)s^ + Ds + 

G(s) • ^ j 

' Is (I -A /'Os + Ds + a 


which is also derived by direct calculations. 


(All) 
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If we further assume that (p(0 takes the form 

<PU) = Y 


where f, is a scalar constant, then by Eq. CA5) 

3 


... 


c = 


35 




= 


and 


6 = — mL (p 


[i^] 


where d) is determined by the normalization condition and becomes in this case 
o • 


Hence, the highest order coefficient of the numerator of Eq. (All) is given 

by 


*«b - l-5(^4)(i) 


Now, we confirm that the poles are not affected by the sensor location, 
while the zeros are significantly affected in view of the equation 


(i - c,(y) 


2 2 
s + Ds + 0 = 0 


with 


a(5 




We observe that as 5 Increases from zero to the maximum attainable value, L, 
a changes continuously from unity to some negative value. Fig. A2 illustrates 
the dependency of the zeros of G(s) on the sensor location parameter a. 
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Chapter 3. Accommodating Model Error in Linear 
Systems Control 

ABSTRACT . This study -is 'motivated by the need for a control design procedure 
which gives to the controlled system three basic properties: a) insensitivity 

to modeling errors (such as from truncated modes and uncertain external 
disturbances) ; b) simplicity of controller (of small order) ; and c) guaranteed 
value of the performance measure of the actual system. This status report 
sketches some progress on the first two goals but does not touch on the third. 
(Actually, property a) is desirable only because it might not be possible to 
mathematically insure property c.) 

Alterations are made to linear regulator and observer theory to accommodate 
modeling errors. The results (some of which are yet unproven) show that a 
"model error vector," which evolves from an "error system," can be added to 
a reduced system model, estimated by an observer, and used by the control law 
to render the system less sensitive to uncertain magnitudes and phase relations 
of truncated modes and external disturbance effects. A procedure is outlined 
to give the observer a "model learning" quality. Two parameters of the error 
system, an "observation window" T, and the dimension of the error system, d, 
are related to parameters of the control problem. By choosing the optimal 
cost as a Liapunov function, the "observation window" (t) of the observer is 
related to the minimum eigenvalue of the closed loop system. Necessary con- 
ditions are given for the solution of x and d. In a rather novel turn of 
events we find that instead of the usual pattern of "given a model, apply the 
control theories" we are using given control theories to help construct a more 
appropriate model of the system. 
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r. INTRODUCriON 


The theory available for analyzing and controlling linear systems is 
quite extensive for LINEAR QUADRATIC problems provided the model Ca set of 
ordinary differential equations) is perfectly accurate. Even if the model is 
in error, one can analyze the effects of these errors post facto . There is 
not, however, a well developed theory which allows one to consider the effects 
of model error during' the control design. 

The exact problem which we would like to solve is the following: Find a 

control policy u°(t) of minimal complexity which guarantees the performance, 
V°. Mathematically, we say it this way: 


THE MINIMAL CONTROLLER PROBLEM : 

Find the control policy, u° = ^(z°,A^ ,B^ ,C^ ,M^ ,t) ,, which though based 

upon the model, , of the physical system, 4 . ' 


x^ - A^ x^ + B^u° , x^ = n, - vector (STATE) 

i 

y* = k - vector (OUTPUT) 


yj _ cJ 
z^ = M^ x^ 


i - vector (MEASUREMENT) 




u = m - vector (CONTROL) 

Physical system, with actual inputs u°, and measure- 
ments z , 

ie Q 

and output quantities , y 


satisfies the performance requirement 


00 

^ 1 1 ly°l Iq I Ctdy^°) 1 ijdt < V° (specified number) 

for minimum n^ , for all expected conditions (disturbance and control 
inputs and parameter variations) . 


In this paper, the word "output" is defined to be the vector of variables 
which we wish to control, (y) . The vector z represents those quantities 
actually measured. 
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Even this mathematical statement of the designer's goal has some vagueness 
that must be Interpreted. The Implication Is that the selection of those 
variables which are to be controlled Cldentlflcatlon of u°) ; selection of 
parameters (A^ , , M^,t); and even the dimension of the model must 

all be determined as part of this "minimal controller problem." Efforts to 
separate the problem Into two separate problems (a "modeling problem" and a 
"control problem") usually yields controllers which either 

(a) suffer from modeling errors, causing performance requirements to 
be violated — 


or, 1 . 

(b) satisfy performance requirements at a cost of high controller 
complexity. i 


The conditions which cause (^a) are sometimes not discovered before "flight." 
The condition (b) often prevents use of modern control theory In an applica- 
tion. Determination of "expected" disturbances and selection of the 
"requirement," V°, can be difficult decisions In the evaluation of a control 
law. 

Presently the "minimal controller problem" can only be solved by trial 
and error. Because we cannot solve this problem and guarantee the cost 
mathematically, we sometimes seek certain precautionary measures by 

1. "worst case" designs (a deterministic approach of us”ing conservative 
conditions In design to add reliability to "flights." This can 
cost us in controller complexity). 

2. Sensitivity approaches (making the system less sensitive to some 
imcertain parameter or disturbances adds some confidence to the 
solution) . 
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3. Stochastic approaches (offer good performance "on the average." 

A particular experiment might violate requirements) . 

In this paper the minimal controller problem (MCP) is not solved, but we 
are guided by the MCP objective to suggest a procedure which allows certain 
parameters and dimensions of the modeling problem (MP) to be related to 
parameters of the control problem (CP). The procedure offers to the con- 
trolled system (in lieu of a MCP solution) 'a certain type of model error 
"forgiveness." To date, however, the cost, , is not guaranteed a priori. 

In section II the modeling problem (MP) is divided into a two phase task: 
a part of the model is determined without regard to the control problem (CP) ; 
a second part of the model is presented in structure but the dimension and 
certain parameters are left to be determined in the CP. Section III presents 
an observer for estimating all the states of the finally selected model. 
Special cases of this observer are discussed which produce the Luenberger 
observer, the Disturbance Absorbing Controller and a Model Learning Observer. 
Finally, in section IV, necessary conditions are given for the selection of 
the control parameters by viewing the optimal cost as a Liapunov function. 

Some concluding remarks appear in section V. 


II. THE MODELING PROBLEM 


Suppose we are given some model of a physical system,, 4 .' Let us label 
this model 


4 


.3 

X 


.3 3 , „3 o 

3 


= A X + B u 

X 

= n_ vector 

^3 3 

3 

J 

= C X 

y 

= k vector 

„3 3 
= M X 

3 

z 

= S, vector 


( 2 . 1 ) 
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For the moment we disregard the means used to obtain We consider 4 ^^ 

* 

to be the actually realized model and we wish to assess the errors which 

accoD5)any 4 ^^ relative to a better model, 

/ .1 .1 1 ^ „1 o . „1 1 

X = Ax +Bu +W, X vector 

J ^ A 1 ^ 

y = C X y = k vector (2.2) 

, z = M X z = Z vector 


f 

The model is /the most general description of which we write ( 

I 

might be used, for example to evaluate performance predictions and controller > 

I 1 

designs which are based upon simpler models). The model ji^ may not haye pro- 
ceeded directly from a prior system description such as However, there 

exists a transformation on jgff such that 4x1^ is a truncation of the trans- 
formed 4fH^. To illustrate this point of view imagine the transformation, 

= T i which has the transformed system description 


^21 ^22 


+ ^ u + 


c\cj.] 


(:;■) 


m'"] / X- 


(2.3) 


T"^B^ = 


t“^a^t = 


A^ A 
A 

^21 ^22 


t"^ = 


, = rC^,Cj, Ml 


(2.4) 


Suppose, for instance, that the parameters actually used in a controller 
design or Implementation are (a3,B^,C^,m3) . 
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Thus X is 


the n^— vector which satisfies 


.3' 

X 


a3 3' 
= A X 


+ 


+ e'Ct) 


(2.5) 


where the vector e' (t) can* be considered a *'model error vector" of the model 
^3 (compare (2.5) and (2,1)). The vector e*(t) evolves from the "error 
system" 


*t = ^ 2 \ ^2^^° + ^2 


e'(t) = A^2*t \ ^ 


( 2 . 6 ) 


which is coupled with and is an integral part of the system description, 
but which is neglected in the system description, Jtdy Equation (2.6) clearly 
focuses the fact that for any model of a physical system there is associated 
with this model an error vector that can be considered as a combination of 
external disturbances and truncated states of a more accurate model. Further- 
more, such an e'(t) exists which compensates for any parameter errors in 
3 3 3 3 

(A ,B ,C ,M ) relative to some intermediate, n^-dimensional, model which the 

designer may have intended to implement. The converse of this statement is 

3 3 3 3 

not true, however. That is, there is no set of parameters (A ,B ,C ,M ) which 
can necessarily compensate for the effects of truncated states and external 
disturbances (e'(t)) which have been neglected in the model. The effects of 
such persistently acting "disturbances" can cause instability in the control 
problem or instability of Luenberger Observers in the observation problem. 

In stochastic descriptions of a system, e'(t) represents correlated disturbances 
which have been ignored. It has been shown by Fitzgerald [1], Price [2J and 
others [3]-[5], that the Kalman filter can diverge when such correlated 
"disturbances" are either neglected or modeled as white. Moreover, even 
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adaptive techniques (which update parameters of the assumed white noise 
statistics) may fall In such cases for reasons stated above. Thus, "structural" 
errors In modeling seem much more critical than so called "parameter errors." 

It is the purpose of this section to propose a model error vector, e(t) , 
which is an approximation to the vector e*(t), and which when added to model 
^ 2 , effectively compensate for modeling errors arising from, 

- truncated states (regarded as "Internal disturbances") 

- external disturbances (which also includes effects of "weak" non- 
linearities) 

- parameter changes 


The vector e(t) is considered to evolve from a much simpler "error system" 
than described by the (n^^ - n^)^^ order system (2.6). The error vector e 
which is added to thusly. 



.3 3 . o . 
Ax + B u + e 

^3 3 
C X 


m3 3 

M X 


(2.7) 


is considered to evolve from the error system 


Y = Dy Y = d -vector 

( 2 . 8 ) 

e = Py P = n^xd matrix 

Now Y(t) is s small dimensional vector which represents a "compression" of 
information in x^,W^ and retains only important effects of large dimensional 
vectors x^(t) and W^(t). Thus, Y(t) is an , artificial vector whose space should 
not be presumed completely controllable because external disturbances are 
not controllable (by definition of the word "external," W^(t) is not causally 
related to x^(t) or u°(t)). As an approximation, we write the error system 
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(2.8) unforced by u . The truncated states may well appear In the measure- 


ments since in (2.3) would not be zero, in general. We therefore allow 
the artificial variables y(t) to Influence the modeled measurement vector, 
z^(t) (via matrix M^) , in the proposed model structure ^2’ 



and the superscript 3 is dropped now that discussion on the hierarchy of 
models is completed.' For our purposes we assume in ^2 (which has the error 
system augmented to model that A,B,C,M are given and we must find the 

error system parameters P,D,M^ which are appropriate for a given problem. The 
reader may well wonder at this point why we wouldn't simply Incorporate the 
knowledge of important disturbance effects in the basic model reduction 
decisions which yield that is a completely adequate model for con- 

f 

trol design. However, we do not usually know what disturbance effects 
(external and internal) are important to keep in the model prior to controller 
design. Performance evaluations tell the control designer, in a trial and 
error fashion , what disturbance effects must fall within the spectrum of con- 
trol authority (i.e. what disturbance effects must be actively controlled). 
Also, parameter uncertainties (whether constant or due to in-flight changes) 
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and variations In the disturbance environment cannot always be reliably pre- 
dicted. This perplexity causes the designer to seek a conservative (l.e. 
"worst case”) design with deterministic procedures or to seek statistical 
procedures which, in essence, can promise good performance only oil the average 
(a particular '^experiment might behave quite poorly). 

It is for these several reasons that e(t) cannot be specified prior to 



Figure 1. Modeling Process. 


*The vector e(t) has the nature of a "correlated" disturbance in the stochastic 
view of the problem and "persistently acting" disturbances in the deter- 
ministic view. 
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Specifying a Structure for P, D 

Since Y(t) is to be selected so that the approximate error system (2.8) 

represents the important effects of the actual error system (2.6), we could 

select the parameters of (2.8) so that 
T 

y*||e'(t) - e(t)||g^^jdt ^ (2.11) 

. ° t . 

is a minimum, where T is the time interval over which good error system infor- 
mation is needed by the controller and g(t) is a weighting scalar. Now, if 

. - _ • -JU. 

we view this intermediate task, as a curve fitting problem in which the vari- 
ables Y^(i) ^re a set of known functions, then the set of coefficients of 

J ' 

those functions which will minimize (2.11) is 


T 

/• 

■ i. 


P = / g(t) e*(t)Y^(t)dt 

•'o 

/ g(t)Y(t)Y(t) dt 
. 0 - 

(2.12) 


One possible disadvantage of least squares methods in modeling problems is 
that the solution which minimizes (2.11), subject to g(t) = 1, can permit ' 
large instantaneous deviatidns between e' and e. This is not a fault of 


least squares theory but can result from poor judgement in deciding what 

/ I I 

to take the least squares of. A wiser choice of the integrand of (2.11) 

I ' 

might include rate terms (e (t) - e(t)) as well. Another option, v^ich we 
pursue here, is to minimize the maximum deviation between e' and e; a result 
which can be obtained simply by choosing an appropriate time weighting g(t) 


and interval T. The result is that (2.12) takes on the form 
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(2.14) 


where a change of variable has been made for convenience 

' t = Y (a + 1) 

gCa) = (1 - 
Y^(a) = cos(i cos ^a) 

The functions V^(o) are called Chebyshev polynomials of the first kind of 
degree 1. To obtain the matrix D, we may differentiate (2.14) and rearrange 

(2.15) 

(2.16) 
dxd 

Now the matrices P and D are determined to within two parameters; the "observa- 
tion window," T, and d, the dimension of the Y vector. A nontrivial task in 
the computation of P Is the determination of e' , the error vector we desire to 
approximate. It may be too difficult to compute e' exactly, by solving (2.6). 
Also the question of what conditions to impose upon to obtain e' might be 
In doubt. Because an observer will be used In section III to continually 
update an estimate of Y(t), we will find that the system will be forgiving of 
certain kinds of errors (uncertainties) in our selection of e'(t). Specific- 
ally, magnitudes and phase relationships of the various modes that comprise 
e* need not be known to guarantee proper convergence of the error system. 


the form to yield 


DY 


where 


D = - 
T 


0 0 0 0 0 
2 0 0 0 0 
0 8 0 0 0 
6 0 12 0 0 
0 16 0 16 0 


*0 
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This notion is made more specific in section IV. For the present, we assume 
that some reasonable facsimile of e* is obtained for use in (2.13). We con- 
clude this section by noting that (2.9) is the form of the "control design" 
model where the parameter matrices (A,B,C,M) are specified at this point and 

specification of the parameters of the error system P(T,d), D(x,d) and M 

Y 

must await solution of the control problem in sections III and IV. 

III. OBSERVER DESIQI 


In this section we are committed to the design of an observer to con- 
tinually estimate x(t) and y(t) , as defined by (2.9). While constructing such 
an observer we will again, as in the modeling problem of the last section, refer 

to a better model to interpret certain observation errors. 

2 

Suppose X (t) is defined as the vector which satisfies model 

I .2 .2 2 . .2 o 

X = n vector 

y = k vector (3.1) 

z = Z vector 


X = A X + B 

2 -2 2 
y = C X 


= 


where 


■ A 

P • 

, = 

■ B ■ 

.0 

D - 


- 0 - 




= [C.O] 


(n+d) -vector 


C3.2) 


2 9 2 2 

Now z = X represents the best combination of x^(t) which corresponds to 

the actual measurement record z°(t). The actual measurement z°(t) does not 

satisfy (3.1), however, because of model errors. In order to write an equality 

2 

relating the state of the model, x (t) , and the actual measurement record, 
z (t) , we must include an error term, e^. 
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’(t) = + e. 


^ o ^2 

= Z - MX 


(3.3) 


where G Is a result of the kind of error e'(t) shown in (2.5). Indeed, if 
z 

is considered to be a model of evaluation quality* then we can write from 
(2.3), 


z (t) = Mx ~ ^ 


(3.4) 


then from (3.3) and (3.4) 

bearing in mind that (3.5) is an approximation and (3.3) is the definition of 

e^. Furthermore, the approximation (3.5) is good only if e'(t) iS small in (2.5) 
3' 3 

SO that X » X . 

We define an (nri-d - £) vector, C(t) , by 


C(t) = r x^ (3.6) 

where the (m+d - A) x (ntfd) matrix F is to be defined momentarily. Now, if 
we augment (3.6) with the actual measurement relation (3.3) we have 




(3.7) 


Assuming is of maximal rank we can choose F so that the Inverse rela- 
tion of (3.7) exists, 

x^ = L.C + L„(z° - e ) (3.8) 

1 z z 


If system tests are not economical, then performance evaluations are usually 
conducted with computer simulations. Here we assume that model is of 
sufficient fidelity to warrant the confidence of performance decisions made 
from it. By this we mean that z^(t) is sufficiently close to z°(t) to write 
z^(t) = z°(t) . 
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where 


»L2 ]. 


r 


-1 


(3.9) 


ff ■ 

Thus, for later reference we have the relations, from (3.9), 


+ L2M^ = I 



-- = I 

M^L2 = I 

(3.10) 

rL2_._ = 0 

= 0 



By differentiating (3.6) and using (3.1), (3.8) we can construct a differential 
equation from which C is considered to evolve. 


C = r a\a + r a\-(z° - e ) + r b^u° O.ii) 

i z z 

Equation (3.11) together with (3.8) form an "observer"; a linear dynamical 

2 

system whose output (relation (3.8)) yield x (if estimation is perfect and 

''2 o 

X otherwise) and whose inputs are the real measurements z (t) and the real 

controls u°(t) , (see (3.11)). Of course e (t) must be available if the esti- 

z 

2 ■ > ■ 

mate of x is to be perfect. Actually equations (3.11), (3.8) offer a struc- 
ture from which a number of different "observers" can be discussed. 

CASE I : The perfect observer 

Equations (3.11), (3.8) and Figure 2 describe a perfect state observer, 

^2 _ 2 2 2 

X = X , if X (t^) and e^(t) are known. (If x (t^) is known, then C(t^) = 
Fx^(t^) and = z° - M^x^ are available, as required,.) Of course, x^(t^) is 
usually not available, else we have no need for the observer at all. 
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(REQUIRES ) 


Figure 2. Perfect Observer. 

(Requires e^lt), x^(tp)). 


CASE 2 : The Luenberger Observer 

If we presume that our model is perfect, then = 0 and (3.11), (3.8) 
reduce to the Luenberger observer for the system 

CASE 3; A Model Learning Observer 

Given the model structure of (3.1), (3.2), the observer could "learn" 

o 2 '^ 2 . 

the correct matrix which would make the "measurement residual," e=z x , 

As 

approach zero. In’ Figure 3, e is used as an error signal to correct M , which 
appears within the gains L^,L2. The physical interpretation of this scheme is 
the following: For a given character of the error system specified by (2.9) , 

this observer will "learn" how the error variable Y(t) propagates to the 


measurements, z =Mx + by finding the matrix which will make the model 


predictions compatible with the measurements, i 


If an algorithm 
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can be found which will converge on the correct then the observer could be 
called a "model learning" observer, in the sense that it learns how the error 
vector e = Py influences the measurements. ' The matrices P, D characterize the 
the error vector e(t) . For Chebyshev polynomial approximation the' matrices 
P, D are given by (2.13), (2.16). 



Figure 3. A Model Learning Observer. 

CASE 4: The Disturbance Absorbing Controller 

If we assume the model is perfect (that is, e(t) E e'(t) in (2.11), such 
that E^(t) E 0) and that the error vector y is totally composed of external 
disturbance effects which cannot be measured (that is E 0) , then the con- 
trol structure of Figure 2 reduces to the so-called Disturbance Absorbing 
Controller whose structure was proposed in (6) and which is developed 
extensively in (7) -(9). 
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One final note is of interest on these observers. Since e is defined 

z 

to be = z° - it is perhaps tempting to construct the approximation to 

A Q 9-^2 

e , namely e = z - M x . If one implements e Ct) this way then the observer 
z z z 

becomes completely disconnected from the measurements z°Ct) , {.see Figure 2) . 

A similar circumstance occurs in Kalman filtering when the plant noise is small 
or zero. Then the Kalman gain can become zero, essentially disconnecting the 
filter from the measurements. Of course, instability is assured in both cases 
since the model which the observer/filter has "learned" contains error and its 
predictions eventually diverge from the real system. Thus, the requirements 
stated in case 1 for the perfect observer can not be met and the attempt as 
outlined in case 1 will guarantee instability because it will not be causally 
related to the measurements. 

The observers described in this section, in conjunction with the model 

» 

given in (3.1), (3.2), can be interpreted as performing the following service 
for the control system: The observer has within it a model of the uncor- 
rected system When measurement records do not correspond to this model 

the difference is curve fitted with a set of Chebyshev polynomials. It should 
be mentioned that the estimates of Y^Ct) are not themselves Chebyshev poly- 
nomials. But over any time interval, (t-T,t), the observer provides y^, (or 
equivalently, Y(t) , since y = e^^y^) , those initial conditions which correspond 
to a best fit of a set of Chebyshev polynomials of degree d over the interval 
(t-T,t), to the actual error vector (MAX[e'-Pe^^y^] is minimized over (t-T,t) 

by the choice of Yq» produced by the observer). For these observers to be 

2 

stable (convergent on x ) the matrix FAL^^ must be a stability matrix and 
(3.10) must be satisfied. The error analysis is incomplete, but the first 
observation is that if C(t) is the (perfect) observer state which produces 
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2 2 '' 

X (t) (via (3.8) and (3.11) with and if C(t) is the state of the 

''2 

physically realizable observer which produces x (t) (via (3.8) with = 0 in 
(3.11) and arbitrary), then the observer error 

A = C(t) - lit) (3.12) 

obeys the differential equation 

1 

A = TA^L^A - rA^L2£^ ^ (3.13) 

It is the purpose of the model learning observer to keep £ small by an 

z 

appropriate choice (perhaps iterative) of M . Noting that e is not zero for 

Y z 

any finite dimensional linear model of a physical dynamical system this is an 
important step in an attempt for better controllers. Luenberger observer and 
Kalman filter techniques assume = 0 (or white) and for that reason stability 
of observer or estimator cannot be assured in general. It is the purpose of F 
to make (3.13) stable and to make (3.7) invertible. This is possible under 
certain observability conditions that will be detailed in further work. 


IV. THE CONTROL PROBLEM 


The model and observer are characterized by the parameters (A,B,C,M,M^,F, 

o ^2 

P(T,d), D(T,d)). The control law to be utilized is the linear form, u = Gx . 
Therefore, the parameters that remain to be determined in this section are 
T, d, r and G, assuming that is determined in section III. 

Selection of T, d: 

Suppose V(x^,t) is a Liapunov function for Equation (3.1), with the 

control 

u° = G x^ (^.1) 

Then, by our previous interpretation of x as an "information window" over 
which good information is needed by the controller, we define X to be a mini- 
mum time constant of the system in the following sense 
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A 

T = MIN 
x^(t) 


VCx^.t) 

-VCx^.t) 


(4.2) 


Now if V(x ,t) is a quadratic form 


V(x^,t) = x^(t)^ ^x^(t) 

• 2 2 T 2 

VCx ,t) =-x^(t)^^x^(t) 


(4.3) 

(4.4) 


T T 

2 '*’ 2 2 ^ 2 

Then we utilize the fact that the ratio = x ^x /x x is bounded above 


and below by 


\cN 




T 

2 2 
X y X 

T 

2^2 
X t/r X 




(4.5) 


to obtain 


T « 




(4.6) 


Since the matrices c/V and 9 are functions of P(r,d) and D(T,d), we could 
find the best T by minimizing (4.6). That is 


MIN 

T 


MIN 

x^(t) 


V 

-V 


(4.7) 


yields T. Rather than solving this problem directly, we will establish an 


upper bound on 


MIN 

x^(t) 


(t) 


and minimize the upper bound. From matrix algebra (10) we can write 

+1/2 




- 1 , 


,,-l| 


^ i 


1/2 


(4.8) 
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where 


^1 


n+d 1 1/2 

Li. 3-1 


1^ 


-li 


n+d 

E "b 

Li,j=l 


2 

ij 


1/2 


(4.9) 


If we choose V(x ,t) to be the optimal cost functional 


2 A 


V(x‘,t) 


aa 

f 


T T 

Cy^ Q + u° R u°)dt 


(4.10) 


and u to be the corresponding optimal control 
u°(t) = ic x^(t) 


(4.11) 


0 = -ic - A" K + K B" R 


2 ^ _ _ 2 - 12 ^ “ 2 ^ 2 


where K satisfies the degenerate Riccati equation 

C“ Q (4.12) 

T T 

then ^ = K, G = -R K, •A' = QC^ + KB^R K, (the observer is con- 

/ 

^2 2 

sidered perfect for the present purpose, x = x ) . Tf we minimize the upper 

' '/ 

bound as specified by (4.8) we have the following necessary conditions for T. 


1? [l 


1/2 




1/2 


3 ||> 


-111 1/2 


9t 




1/2 




(4.13) 


.-1 


Equation (4.13) is quite difficult to solve and we must show that JV exists. 
1/2 ••I 1 / 2 

To obtain 1 | ^ | 1 , | \<A | | in a form explicit in T requires a Kronecker 

product (11) solution of 


-L(T,d)D(T,d) + [KBR V - A*^] L(T,d) = K P(x,d) 


(4.14) 
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for L and of 


-F(T,d)D(T,d) - D^(T,d)F = + P^(T,d)L + L^P(T,d) (4.15) 


for F, where L and F are partitioned parts of K in (4.12), 


L(T,d) 


L'^(T,d) F(T,d) 


(4.16) 


Therefore 


E lU. 


(4.17) 


where £. are columns of the matrix 
\ 1 


[;]- 


1*^2’* “^d^ 


(4.18) 


' Further work is required to determine If T can be solved from (4.13) and 
whether It Is unique. 

Assuming that the norm of ^ Is a convex function of d, we will select d 
" to satisfy 


T 

A I 
d d 


a E 


i !■ . 


(4.19) 


for small a » .01. This In effect forces a minimum degree of observability of 
"error system" modes (l.e. so that the d^^ (and final) mode of the error system 
Is say, 1% as observable In the output as - all other "modes of the error system) . 
Mathematical justification of this Interpretation will he offered later. 
Choosing a to be as large as, say .01, serves to constrain the order of the 
error system. Otherwise d might- be quite large (to obtain arbitrarily small 


errors In fitting e(t) to e'(t)). 
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V. CONCLUSIONS 


An "error system" is augmented to a system model which is intended for 
controller design. The parameters of the error system are determined in such 
a way that the observer, which is used to estimate the. states of the model and 

error system, serves to fit actual model errors (the difference between real 

■ - • 

measurements and those predicted by the model) with Chebyshev polynomials so 
that the maximum modeling error is minimized over an "observation window," T. 
In this view of the operation of the observer, the system can perform, in 
essence, "adaptive curve fitting" of internal and external disturbances with- 
out recourse to "adaptive" techniques. Alternately, in the mode of operation 
in which the observer is adaptive, the observer could be called a "model 
learning observer" in the sense' that it learns how the model errors Influence 

the measurements. To accomplish this an algorithm must be found for changing 

A ■ o* 2 

a certain matrix, M , so that the measurement residual, e = z ' - z , Is 
Y . ^ ' 

driven toward zero. This adaptive feature has not been completed, however. 
Certain parameters of the error system model must be determined simul- 
taneous with control policy design. Specifically, the "observation window," 

T, associated with the error system is shown to be related to the minimtnn 
eigenvalue of the closed loop system. This relation is accomplished by 
viewing the optimal cost as a Liapunov function. The dimension of the error 
system, d, is determined fr 9 m a required degree observability of the modes of 

the error system, although this feature has not been completed and only heces- 

< ■ 

sary conditions are shown for T and d. 

■ As a speci£il case of the observer derived herein, the Luenberger observer 
is Obtained. As another special case the method reduces to the Disturbance 
Absorbing Controller reported in [9], [8] and [7], The present method has 



several advantages over the Disturbance Absorbing Controllers. Most signifi- 
cantly, the effects of "Internal disturbance's" (truncated’ modes ) are not 
Ignored (l.e. they are not treated simply as external, unmeasurable 
disturbances) . . c 
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Chapter 4. Hybrid Coordinate Method Using 
Assumed Mode Shapes for Elastic Contlnua 

ABSTRACT ; The hybrid coordinate method provides equations of motion of 
minimum dimension for a spacecraft with flexible appendages. Instead of 
the usual finite element approach. In which mode shapes are calculated from 
equations of vibration of the finite element assembly, this chapter provides 
an alternative formulation using assumed mode shapes. This proves useful 
for a class of simply modeled appendages for which mode shapes are provided 
by an outside agency, or are otherwise known. The results are shown to be 
compatible with the finite element formulation, as previously described. 
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I. INTRODUCTION 


The hybrid-coordinate method provides equations of motion of a space- 
craft with elastic (flexible) appendages. The appendages are modeled as an 
interconnected set of small rigid bodies interconnected by massless or mass- 
ive elastic bodies (finite elements). From the Newton-Euler approach, the 
equations of motion for each finite element and the rigid body portion of the 
spacecraft are formulated. The Introduction of an appropriate coordinate 
transformation allows the finite element equations to be represented as 
decoupled vibration equations, which Involve mode shapes and modal coordinates. 
Since the vibration equations have been decoupled from each other, significant 
truncation of the higher order mode shapes can be accomplished. This leads 
to a set of equations where rotation of the rigid body portion of the space- 
craft is coupled to the vibration of the flexible appendages. These equations 
are of great practical use because the truncation procedure has significantly 
reduced the number of degrees of freedom of the system without substantially, 
sacrificing the fidelity of the results. 

the purpose of this chapter Is to provide an alternative formulation 
for the hybrid coordinate method using assumed mode shapes. This approach 
will prove useful for simply modeled appendages. For these the mode shapes 
can be determined from a continuum analysis using partial differential equa- 
tion methods. The truncation procedure Is accomplished at the outset by 
erimlnating th^hlgfi« order modes of vibration. The equations of motion 
are formulated using a Lagranglan approach and the coordinate transformation 
Is accomplished using the assumed mode shapes. The resulting equations of 
motion are then seen to be compatible with those arising from the finite 
element method. 


/:/// ) J ; 
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- Figure 1. System Diagram. 

Center of mass of undeformed system (body fixed) 
Position of 0 at rest (inertially fixed) 

Center of mass of rigid body 
Center of mass of undeformed appendage 
Total System mass 
Appendage mass 

\ 

Appendage mass/length 
Connection point of 'appendage 





II. MODEL 


The following derivation of equations of motion. uses a model comprised 
of a central rigid body and a flexible cantilevered beam. Extensions to 
'several appendages of arbitrary configurations may be made from the results 

of this simple model. The undeformed position of. the appendage Is taken to be 

constant relative to the rigid body. The transformation between the two Is 
Included In the derivation to facilitate the extension of the equations to 
cover a driven appendage. The angular rotations are assumed to be small as 

are the translational displacements. A diagram of the model Is shown In Fig- 

ure 1. To summarize, the assumptions used In the following derivation are: 
s rigid body with cantilevered beam, 

• beam rest position constant relative to base, 

• small translations and rotations^ 

No orthogonality requirements have been placed on the assumed mode . 
shapes. The vibration equations are therefore coupled. Further coordlnatd\ 

■ . i ' , ■■ 

transformations may be employed to decouple the vibration equations or to . 
achieve vehicle normal modes,, but the truncation procedure does not require 

I 

this as It does with the finite element procedure. 

j 

The vector bases emplojyed In the derivation are: 

1 • 

{i} : Inertlalljij fixed basis 

• I 

{b} : Basis fixed In the rigid body 

{a}. : Basis fixed in appendage prior to deformation 


where 


{b}' =. [9]{i} 

[01j= (E - 0) 

I 

la>: = [c]{b> 

I 

[c]j = constant for an undrlven appendage. 



’ 0 -0- 

®2* 

for small rotations 0 = 

I \ 


. 1 ^ 
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For the undeformed system, the location of the center of mass is 
defined by 

/ 

p dm = 0 


/h 


SYS 

where p Is the generic position vector from the center of mass to the dif- 
ferential mass element. Evaluating this expression leads to 

-(v<r-M) L + m(r + r^) = 0 

where the quantities are shown In Figure 1 with L being the location of the 
rigid body center of mass and r^ the appendage center of mass. 

The dyadic of the undeformed system Is defined by 

D*-Q“ ^D° 

RB ^AP-u 


n°=n® + -M) (L*L U-LL) 

'-'rb ^RB 

' 0° ■ * " [(—*)■ (— *> - ■ (S*Ja)(£*Sa)] 

AP-u AP-u , •' 

III. THE LAGRANGIAN OF THE SYSTEM 
The kinetic energy of the system Is 

•|y* V • V dm + jJ" V • V dm 


T 

SYS 


RB 


APP 


where V Is the Inertial velocity of a generic mass element. The kinetic 
/ . - . . 
energy for the rigid body yields • - , - 


1 y* Y • V dm = Y + Y w * D * y 

■ RB 


RB 


where tJif Is the system mass, M the appendage mass,, and B Is the center of mass 
. . ■ . .. .. 
of the rigid body with 
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• • 

Rg ® Z “ cu ^ L« 

Expanding this expression and switching the reference point o^ the Inertia 
dyadic to the system mass center gives 


iy* V • V dm - ) Z«Z+ -i o) * • (a»xt)j , 


RB 


The kinetic energy of the appendage Is 


where 


ly V . Vd..iyR^ . dm 
APP 


R^“Z + R + r + u 


o o 


• . . u u / > 

^«Z + u + R + r + wx / R+r \ . 


Vector differentiation with respect to the rotating reference frame is denoted 

0 o' ' 

by the "circle" above the vector. Here R and r are zero since they are fixed 
in the frame. The "dot." denotes differentiation relative to an Inertial refer- 
ence frame. Expanding the expression and making use of the dyadic of the unde- 
formed appendage about the system mass center (fl Wields. 

'^AP^u' j 

hfl' “ I f y • • y + 7 y*y • ' 

+'Z • y* u dm + Z • u) X ^IR + r^^j 


APP 


APP 


/[^ • y X (R+rjJ dm 


APP 


Combining the terms for kinetic energy and eliminating terms produces 
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T 

SYS 


J 40 • [] *40 


+ ^ J“ u»u dm + Z • j' u dm 
APP APP 

w*Rx^ u dm + u • J' (rxu) dm. 


\ 


APP 


APP 


The inertia dyadic of the tindeformed system about the system center of mass Is 
Jt 


D' • J 


Tbe' center of mass expression eliminated the term In the kinetic energy 
/containing Z and oo. 




From beam theory, the strain energy of the appendage Is 


- \ , Up “&)■&) 

I I I 

The Lagranglan for the system Is then 


dr 




h-2'^ 

Z • J' u dm + (j) • R J' u 


APP 

dm 


APP 


APP 


/ 2 \ / 2 \ 

Di • f (rxu) dm - ^ r El [ ^ ] dr. 
J ~ ~ •' \8r / \3r / 

Ann ' r \ f 


APP 


APP 


/ 


The formation of the above has ^sumed that the undeformed appendage Is 
fixed relative to the base • The next step to be taken Is to 

assume small angle rotations and represent the Lagranglan In matrix form. The 
following matrices are used: j 


Z = {1}^ {Z} 


U) - {b}^ {0} 


(small rotations) 
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u = {a}^ {u} 

Q*- {b}^ I* {b} 
R = {b}'^ {r} 
r = {b}^ {r} 



Retaining second order terms in the Lagrangian produces 

^ {z} + i {0}'^ I* {0) / 

+ y* {u}^ {u} dm + {z}[C] ^ {u} dm! I 
APP APP 

+ {0}^ ^ dm + c y* r {uj dm 

APP APP 

1 /* I'l Ti It ! 

- / Ei {u } {u } dr. 

APP , 

I 

I 

153 

I 



Distributed coordinates are Introduced by the coordinate transformation 

n 

u(r,t) = 53 $^(r) (t) 

" 1=1 

where n is the number of modes used to represent the displacement. In matrix 

form, the transformation Is 
n 

{u} = 2 

1=1 

” I I ••• I 

- [4>Un} 

where [ip] Is a 3 ^ n matrix with each column corresponding to a mode shape 
and {r]} contains n modal coordinates. 

This coordinate transformation yields the Lagranglan , 

5? = Y ^ {Z}’^ {Z} + -| {0}^ I* {0} 

+ 1 {fl}^ X 2 {fi} + {Z}^ C {fl} 

+ {0}^ C(R Xj^+X 3 > {p} - j {p}^ X^ {p} 

where 

X 3 = A [<|)] dm 
APP 

X 2 = /* [<|)] dm 

APP 

X 3 = ^ “r [p] dm 
APP 

/ II itT fl 

El [p [p ] dr 

APP 
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The matrices and are of dimension 3 x n while the matrices X2 and iX^ 
are symmetric and of dimension nxn. The Lagrangian depends on' n+6 gj'eneraiizdd 
coordinates. Six coordinates describe the translation and rotation lof the 

... - ' I 

undeformed system and n modal coordinates describe the displacement (from irest 

' I ! i 

of the flexible appendage relative to the rigid base. 


I -J 


IV. EQUATIONS OF MOTION 


The equations of motion for the system may now be derived from' the 
Leigranglan in the traditional manner. The resulting n+6 equations jmay be 
represented in mat^x form as 

.^{Z} + Xj^ {ri} = 0 ! i I 

I* { 0 } + (r X^+X3 ){*n} = {T} , ' 

X2 {ri} + {n} “U-X^ {z} + (x^ R-X3){e} 

where^jCri-'ls'^^ externally applied torque. The first matrix equation may j be 
used to eliminate the translation from the vibration equations. This pro-| 

duces n + 3 equations of the form | 

■ ■ ! j ■ 

I* { 0 } + (r X3 + X3){n} = {T} ' ' 

(x2 X3 xj{f 1 } + X^ {n> “ (x^ R - x^) {g} 

' ■ I 

The matrices that provide coupling between the rotation and vibration 


in each equation may be seen to be transposes of each other. The 
can be written as 

, I* {0} ^ 6”^ {ri} => (t) 


equations 


(^2 ^1 ^ 1 ) + X 4 { n ) = 6 { 0 } 


where 


X3 R - X3. 
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V. COMPATIBILITY WITH FINITE 
ELEMENT EQUATIONS OF MOTION 

The equations ot motion derived from the continuum analysis are similar 
In structure to those derived from the finite element analysis (Ref, 1), The 
differences between the two appear In the assumptions made in the continuum 
analysis t 


• No orthogonality properties 

• No differential rotation of appendage mass 
elements due to deformation, 

... 

Orthogonality properties can be applied to the contlnuiun analysis 
vibration equations by a suitable coordinate transformation. The orthogonality 
properties are not needed to permit truncation as Is the case In the finite 
element analysis. 


The equations of motion from a finite element analysis are /shown by 
equations (287) to (289) of Reference 1. 

I* 0 - 5'^ fi = T 

n + 26ofi/+ a^n = 60 

/ 




The overbar Indicates truncation. If the damping Is eliminated and 
the orthogonality condition relaxed (after truncation) , the equations become 

I 0 - 6^ n = T 


where 


i m' (j) n + k' <J) n = t 0 

m’ .M(E-EEo'2Eo’'»'^«'j 


Stlffnes- matrix. 
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Normally, the coordinate transformation ^ Includes mode shapes with 
translation and rotation of the finite elements. To agree with the continuum 
analysis, no rotations of the finite elements will be allowed. The coordinate 
transformation ^ will then be a 6n x N matrix represented by 


where 


$ = [{4.^} {4>^} ...{<!)“}] 


r il 

{^Dl} 




(♦J } 

n 

• 0 • 


With the above limitations, the matrix multiplication can be performed 
in the finite element equations and the terms may be compared with those 
from the continuum analysis. 

For the augmented mass matrix, the finite element analysis results in 

IT 


n 


($V) = "V 

' ij k-1 


4> 


4.J 


n n . 

Z Z 

A=1 k=l \ 'JT 


<P 


rr 


This. is compatible with the result from the continuum analysis 

,T 

m(r) j 4>^ 

'ij APP 


(h -hr *1) '/ 

'A A w U J 


■// 


dr dr. 
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For the 6 matrix, the finite element analysis results In 


[ 5 ], 


1=1 “i 


i . 




/w . 

(R . r,) 


This is c^patlble with the result from the continuum analysis 
[5]j =y* m(r) ^R + r jdf ^ 

' • I 

‘ . 1 

■ i ' 

Thus,! if t^e number of finite elements were Increased without limit, the 

I ! 

f inlt'e element equations would be identical to the continuum analysis 

i ' ■ 

equatlonsi. 

■ ; ■ ■ I ' ■ . 

: ' VI. CONCLUSION 


With the foregoing results it becomes possible to accomplish a hybrid 
coordinate dynamic analysis for a system with appendages defined only In 
terms of modal data based on a continuum analysis. This is a necessary step 
before we can simulate an LST vehicle with FRUSCA solar panels defined by the 
modal data provided by HAC. 
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